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SUMMARY 


This report describes the last three years work of a five-year 
NASA -sup ported study of the interacting effects of gravity and geo- 
metry on the peak minimum and pool boiling heat fluxes. The work 
has generally sought to clarify this interaction by seeking an under- 
standing of the hydrodynamic mechanisms which dictate it. Accordingly, 
this report emphasizes that portion of the last three year's work which ^ 
has extended the hydrodynamic theory. 

The work in this area, embracing two doctoral dissertations, two’ 
masters theses, two special undergraduate projects, and a variety of 
other studies, represents a major extension of the theory as a whole. 
The major elements of the report include: 

1. ) An extension of the hydrodynamic theory of the peak heat 
flux for flat plates of both finite and infinite size, supported by a 
great deal of new data. 

2. ) A general theory of the peak heat flux on finite bodies of 
various configurations. The theory is applied to a variety of particular 
configurations and supported, in these cases, by much new data. 

3. ) Viscous theories of the Helmholtz and Taylor stability 
problems. These theories are born out in great detail with new data 
for film boiling and for the peak heat flux, in viscous liquids. A 
viscous theory of the peak heat flux is developed and verified experi- 
mentally. 


4. ) A hydrodynamic theory, fully verified with new data, for 
high current-density electrolysis. 

5 . ) An examination of the deterioration of the conventional 
boiling curve, that takes place when size or gravity are reduced to 
the point at which inertia ceases to be of importance. The result is 
a new phenomenon for which a new set of analyses are developed and 
corroborated with much original data. 
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I . INTRODUCTION 


During the past five years, the NASA Lewis Research Center has 
supported an ongoing study of the hydrodynamically dictated peak and 
minimum pool boiling heat fluxes at this laboratory. The results of this 
effort have ranged broadly over many related areas and they are reported 
in 28 project-related publications [1 to 28] 1 . 

The general functional objective of the study was to learn how 
gravity influenced the peak and minimum pool boiling heat fluxes in var-» 
ious configurations. The primary contribution of the work has been to 
broaden the hydrodynamic mechanisms of these heat fluxes, as proposed 
by Zuber [29] and Kutateladze [30], into a general body of theory which 
can describe a large variety of situations. These mechanisms are the 
combined result of gravity, inertia, and capillary forces. 

The initial stages of this program were summarized in the NASA 
Contractor Report [7] of the first two year's activity (Sept. 1967 to 
Sept. 1969). The first extension of the grant (Sept. 1969 to Dec. 1970) 
was reported by a brief letter listing the publications of that period. 
Therefore, this final report will detail the work of the first and second 
(Jan. 1971 to May 1972) extension periods as they relate to the hydrody- 
namic theory. Those efforts of the project which do not bear directly 
on the hydrodynamic theory will only be summarized briefly at the end 
of this chapter. 

Hydrodynamic Theory 

* 

The present report presents the following extensions of the hydro- 
dynamic theory of the peak and minimum pool boiling heat fluxes, which t 
have been made since the last contractor report [7] was presented: 

Chapter II describes the hydrodynamic theory of the peak and min- 
imum pool boiling heat fluxes on infinite flat plates. The theory, as 
originally formulated by Zuber, is altered in the light of recent experi- 
mental evidence and tested against new data. Chapter III is devoted to 
the adaptation of the theory to a variety of finite heaters. Many new 
experimental data are offered in verification of the resulting predictions. 

The variable of viscosity is introduced in the IVth Chapter. Here 
the Taylor and Helmholtz stability theories, which underlie the hydro- 
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dynamic processes, are formulated with consideration of the role of such 
variables as liquid depth, vapor density, and finally liquid viscosity. Ex- 
perimental data are provided to verify the theory. In Chapter V, these 
results are used as the basis for predicting the peak heat flux in viscous 
liquids. Original q ax data are given in verification of the prediction. 

Chapter VI deals with the application of the hydrodynamic stability 
concept to high current-density electrolysis processes. Experiments and 
theoretical considerations show that the peak electrolysis gas flux does 
not correspond with the Helmholtz-unstable gas removal flux, as is the 
case in boiling. Rather it occurs at the Moissis-Berenson point of first 
hydrodynamic transition. A region of film electrolysis is also identified. 

An important limitation on the hydrodynamic theory is its failure in 
situations where surface tension completely overbalances inertial forces. 
This occurs when heaters are very small or gravity approaches zero. 
Chapter VII treats the boiling process on very small wires in an attempt 
to learn how the transition to low-gravity behavior occurs. The peak 
and minimum heat fluxes are found to vanish completely. They are re- 
placed by a heat flux vs. temperature curve that rises monotonically 
through a mixed natural convection and film boiling regime. 

The objective of learning the role of gravity has been satisfied in 
that every theoretical result given in Chapters II through VII incorperates 
gravity explicitly, and every experimental correlation incorperates the 
variable of gravity in the nondimensionalization. 

Other Results 


« 

A great deal of the work done on the project has been essentially 
omitted from the present account since it would add unduly to the length 
of the report. We shall therefore describe the deleted material briefly 
before moving on to the hydrodynamic theory: 

Condensation: The practical problem of refluxing boiled liquids within 

the centrifuge led to theoretical work on condensation and natural convec- 
tion in variable gravity fields. A first paper on laminar condensation 
[10] showed how to rewrite Nusselt's equation in a form that will apply 
when gravity varies along the condenser, and which will handle a variable 
radius of curvature as well. The equation differs from Nusselt's equa- 
tion only in that it replaces gravity with an "effective gravity" which 
takes account of the arbitrarily variable gravity and radius. 
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A second paper on laminar condensation [23] showed how an ar- 
bitrary variation of wall suction could also be incorporated into this sort 
of analysis . 

Convection: A second problem that was motivated by the reflux con- 

densor design was that of predicting natural convection heat transfer in 
a variable gravity field. This led us to formulate a general integral 
method for predicting natural convection in a variable gravity field [22], 
The method was applied to a variety of rotating, and curved -surface, 
heaters. A series solution was developed for the full differential equ- 
ation to provide verification of the accuracy of the integral method. 

Fin design: A brief note was written to show how to design a "per- 

fectly effective" fin [8]. By boiling a coolant on a fairly short, iso- 
thermal fin, near the peak heat flux, an extraordinarily high heat re — 
moval rate can be achieved. Design relations are given for such a fin. 

Electrolysis bubbles: An important task within the last grant from 

NASA was that of exploring the analogy between gas formation and remo- 
val in electrolysis and boiling. The task was divided into two portions: 
a) a study of the analogy between electrolysis and boiling bubble growth, 
and b) the study of the analogy between the hydrodynamic stability mech- 
anisms in electrolysis and boiling which is treated in Chapter VI. 

The bubble growth analogy is dealt with in references [11] and [17]. 
In these studies, the equations for bubble growth are non-dimensionalized 
in such a way as to make them identical for either electrolysis or boil- 
ing. Comparisons are then made between predictions of electrolysis 
bubble growth and measurements of vapor bubble growth. 

Computer coding : A very large number of data had to be cast into 

various schemes of non-dimensionalization. For convenience of handling 
the data, a computer code was developed. The complete data of the 
project were put on cards and the equations for the physical properties 
of the boiled liquids were assembled in a subroutine. The routine makes 
it possible to call forth plots of any appropriate data in any desirable 
scheme of non-dimensionalization. This code is described fully in refer- 
ence [13], 
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II. PEAK AND MINIMUM HEAT FLUXES FOR 
INVISCID LIQUIDS ON INFINITE FLAT PLATES 


The Original Hydrodynamic Theory 

The hydrodynamic theory of Zuber showed rationally in 1958 [29] 
why the older correlative equation of Kutateladze [30] for the peak heat 
flux, q max , on a horizontal flat plate: 

1/2 1/4 

q = (tt/24) o h [ag(o f -p )] = q O) 

max^ 9 fg '9 max Z 

2 

was correct . He formalized Kutateladze T s suggestion that the peak 
occurred when the vapor outflow from the plate became enough to obstruct 
the returning liquid flow: Figure 1 shows his physical idealization of 

the process. 



Fig. 1 Zuber's vapor jet configuration for boiling from a 
horizontal flat plate heater 


The square array of jets is spaced on the "most susceptible" or 
most rapidly growing Taylor unstable wavelength. This is sometimes 


Symbols not explained in context are ones in common use. They are 
explained in the NOMENCLATURE section. 
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called the "most dangerous" wavelength; hence the symbol, X , is used 
to identify it. The reason for this choice is that the vapor generated at 
the plate is lighter than the liquid above it and the configuration is 
inherently Taylor-unstable. 

Zuber had access to no experimental data of a kind that would tell 
him whether the interface would collapse into jets on the most susceptible 
wavelength, the "critical" (or shortest unstable) wavelength, \ , or some 
other wavelength entirely. Hence he was forced to guess that°the col- 
lapse would occur somewhere between X and x^ = /3 \ . We shall pro- 
vide a good deal of evidence related to £his point in Chapter IV. 

The theory of Bellman and Pennington [31] gives X for a hori- 
zontal inviscid interface between a liquid (f) and a gas (g^ 1 with surface 
tension, a, as: 


= /*3 2rr 


V C /9 (p f “ 


P ) 

g - 


= /3 X 


( 2 ) 


where the subscript, F, is added to remind us that this result is restric- 
ted to an undisturbed interface that is flat. 


The jet pattern is now known to be well-established at heat fluxes 
far below d max * In 1962, Moissis and Berenson [32] showed that the jets 
came into being as the result of a hydrodynamic transition from isolated 
bubble behavior to continuous columns. They gave an expression for this 
transition and corroborated it with experimental data. Thus Zuber's pro- 
posal that the vapor generated at the wall is carried away in these jets 
was shown to be true over much of the range of nucleate boiling. 

The jets continue to serve as the vapor escape route until the vapor 
velocity is high enough to make them Helmholtz unstable. This instability 
causes the entire process to collapse and it defines the maximum attainablt 
heat flux in nucleate boiling. We can then write the peak heat flux as 
the maximum latent energy transport from the surface: 


q = h (A /A ) U 
max g f g v j hr H 


(3) 


where A. and A are the cross-sectional area of the jets and the area of 
the heater, and U is the vapor velocity for which the jets become 
Helmholtz unstable. 


Equation (3) is perfectly general. It applies to any configuration 

if the right values of A. /A. and U are used. To evaluate A /A for a 

J h H j h 
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flat plate one assumption must be made, namely that the radius of a jet 
is > /4 where K is the actual wavelength — be it x^» X^, or some other 
value. This assumption has been supported by some experimental evid- 
ence (see e.g. [33] or information presented later in this chapter). 

Using it we obtain: 


A /A 
j 



it/ 16 


( 4 ) 


The critical velocity, U , can be shown (see e.g . [34], pg . 462 
or [12] to be 


U 


H 


/ 2- c/o 


/v i i 

g h 


( 5 ) 


where \ is the disturbance which amplifies and becomes Helmholtz- 
unstable^ This result, like all of the results we shall present in this 
chapter, is restricted to system pressures that are not near the critical 
pressure. This means that A is « 1 and that the liquid inflow 
velocity is negligible with respect to the vapor outflow velocity. 

Finally, we must determine the appropriate value of \ . Zuber 

guessed the critical Rayleigh wavelength to be the dominant existing dis- 
turbance in the jet. A vapor jet in a liquid collapses in capillary waves 
just as a free liquid jet in a gas does. Rayleigh (see [34], pg . 473) 
showed that this wavelength is equal to the circumference of the jet, or 
r, X/2 in the present case. 

Finally, Zuber substituted equations (4) and (5) in equation (3). 

Then, using equation (2) for X in X = ~X/2, he obtained equation (1). 

In formulating expressions for q and q , Zuber actually be- 

gan by thinking about the transitional Scnfing regTrne and he conceived 
the peak and minimum heat flux transitions as limiting points of this 
regime. Throughout the transitional boiling regime, vapor blankets the 
heater in greater or less degree. The liquid spills through the vapor 
from time to time, either touching the surface or at least drawing very 
close to it. As the surface temperature of the heater is reduced, more 
liquid contact can be sustained and the heat flux is improved. Finally, 
the Helmholtz-stable jet is re-established at a heat flux equal to Q ax * 

An important feature of the transitional boiling regime is that it 
involves a vapor blanket beneath the liquid bulk. Thus this blanket must 
constantly collapse in Taylor-unstable waves to release vapor and admit 
liquid. The vapor escape paths will then be arranged on a rectangular 
grid of Taylor-wavelength spacing everywhere in the transitional regime. 
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The minimum heat flux is reached when the heat flux matches the 
minimum rate of latent heat removal that will sustain the regular form- 
ation and collapse of Taylor waves in the interface over the plate. If 
the heat flux is anything less than this, the interface will collapse onto 
the plate cooling it and re-establishing nucleate boiling. The minimum 
appears to occur in practice after the Leidenfrost point (beyond which 
the surface is too hot to permit direct liquid contact) has been reached. 

The minimum heat flux, q . , can then be written as 

rmn 


, rvolume of~ ."bubbles per cycle - '; [minimum wave] 
H min K g fg La bubble J Ifor unit area J [frequency 

Figure 2 shows the vapor removal configuration as it typically appears 
in reality. 



Fig. 2 The film boiling of acetone near the minimum heat flux. 
0.646 mm dia. wire. 


Zuber assumed the radius of departing bubbles to be \/4. This, indeed, 
was the motivation for his assumption that the jet radius is X/4 and it 
has subsequently been verified experimentally [3CT]. The number of ^ 
bubbles released per cycle in a unit area is 4/),^ for the flat plate. 


Zuber erred on this point and wrote 2/\ (see [12], pg . 25). However, 
his number turns out to be correct for adegitimate 3-dimensional interface. 
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Finally, the minimum wave frequency poses the most serious dif- 
ficulties since most of the growth of the wave up to the point of bubble 
release is non-linear. Zuber attempted to guess the relation of the 
average growth rate to the initial linear growth rate of the wave. 
Berenson [35] subsequently objected, pointing out that there is insuffi- 
cient information to predict this average. Additional data is required 
to get the minimum frequency. 


Berenson argued that the minimum wave frequency is directly pro- 
where 


portional to -iju , 
^ max 

of the linear wave. 


uj is the frequency of exponential growth 
and a pure imaginary number. Using these results 


and writing the volume of the bubble as (4rr/3)X ( _|/64, we get 


^ . = constant p h X , f iuo ) 

nurij. f fg d^ max' 


< 7 ) 


The theory of Bellman and Pennington [31] includes the expression for 

(iuu ) 
max F 


(-U 


,)^= I 


/2tt rgrV 


. 2 

4tt o 


1/2 


max'F ; X , l (p + p ) 2 

d VH f V X d (p f + C g ) 


combining equations (2), (7), and (8) we obtain at last 


( 8 ) 


,cgCp f -0 

q = constant c h / — - 

min g fg 2 

F (P*+P„) 


( 9 ) 


The constant in equation (9) must really be determined experi- 
mentally, Zuber originally suggested a value of 0.177 which is too high. 
Berenson’s experiments subsequently fixed it at 0.09 and that is our best 
information to date. 


While we have gone into detail in describing Zuber^ enunciation of 
the hydrodynamic theory, we shall, for completeness, include a brief 
chronological review of the overall development of the theory prior to the 
period covered by this report. 

1948: Kutateladze [30] proposes the choking mechanism of the peak 

heat flux, and uses dimensional analysis to obtain an empir — 
ical expression for q on a horizontal flat plate heater. 
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1956: Borishanski [36] extends the correlation of Kutateladze to in- 

clude the effect of viscosity. For large flat heaters he 
finds that 


max 


. 2/5-' 

= 1 + 30 . 5/N q 


F viscous 


max r 


( 10 ) 


where N, the "Borishanski number” is: 


N 




1/2 



(ii) 


1957: Chang [37] observes that the Taylor instability is basic to 

natural convection from heated horizontal plates, and to film 
boiling as well. 


1958: Zuber and Tribus [29], [38] propose the analytical model of 

the burnout process that we have discussed here. 


1962: Berenson [35] reconsiders Zuber’s theory both experimentally 

and analytically. He supports the Q m ax theory with an experi- 
mental verification of the predicted insensitivity of q max to 
surface condition. He also shows, both analytically and expe- 
rimentally, that the constant in equation (9) must be determi- 
ned empirically and that it can be altered by surface effects. 

1963: Lienhard and Wong [33] show how to adapt Zuber’s q^jp, pre- 

diction to a horizontal cylinder configuration. The work in- 
volves a prediction of the radius-dependent Taylor-wavelength 
for this configuration. In 1969 [5], Lienhard and Sun provide 
an extension of these ideas and a discussion of some of the 
limitations on them. 


1964: Bobrovich, Gogonin, and Kutataladze [39] use dimensional 

analysis based on the hydrodynamic theory to show that q ma x 
data for complex geometries can be correlated with: 


q / a 
max m ax_ 


= fn (L f ) 


( 12 ) 


where L T is a scale parameter based on a characteristic 
length, L: 
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L' 


L V 9 CPf - P g )/cr 


(13) 


This work was also done independently by Lienhard and 
Watanabe [40] in 1965, who also provided some theoretical 
justification for equation (12). 

1970: Lienhard and Keeling [4] discuss the inclusion of the Bo- 

rishanski number in correlations in general. They 

show that when heaters are arranged so that the rising 
vapor induces a secondary flow of liquid into the jets and 
columns, the influence of N can be very great. 

1970: Sun and Lienhard [6] make a successful hydrodynamic pre- 

diction of q^ax for horizontal cylindrical heaters. 


Modifications of the Theory 


Zuber's theory for involved assumptions which can be ques- 


tioned in the light of things that have been learned during the intervening 
years. These questions will lead to a modified prediction of the peak heat 
flux on a flat plate, which we designate as q as opposed to Zuber's 

m 3.X r— - 


value. 


q 

m ax_^ 


Since experimental data [33] verified Zuber's assumption that the jet 
radius, Rj, equals X^/4 for horizontal cylindrical heaters, we shall accept 
this assumption as being true for infinite flat plates as well. 


The length of the Helmholtz unstable wavelength in the vapor jets is 
very hard to know a priori. It was found in the study of q max on horizon- 
tal cylinders [6] that Zuber's use of the Rayleigh wavelength is probably 
valid as long as the heater (and, with it, the jet) is small. Indeed the 
Rayleigh wavelength can be seen clearly in the jet leaving a small sphere 
[15]. But on larger bodies, from which the vapor leaves in a larger jet, 
the Rayleigh wave becomes much longer than X^. And ^ is the dominant dis- 
turbance in the liquid-vapor interface. This wavelength is picked up from 
the interface and carried into the jets in a rather advanced state of develop- 
ment. 


Accordingly the jets may collapse at a lower vapor velocity [ cf . equa- 
tion (5)] due to Rayleigh instability but only after the vapor has been carried 
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far out of harm's way. The collapse will not lead to burnout. At a 
somewhat higher vapor velocity, however, the shorter disturbances 
of length, , will be triggered and the process will collapse at the 
surface. F Thus we pick X^ = instead of 2nR., for use in 

equation (5), when we deal with the relatively large jets above a flat 
plate. Using equation (2) for X d this gives 


U 


H 


:F 


i( p f -p )/3 Pg ‘ 


(14) 


The area ratio A^/A^ will be tt/ 16 regardless of the exact value of X. 
Substituting equations (1), (4) and (14) in equation (3) we then obtain 

q 

max,_ 

— = 1.14 ( 15) 

H max z 


Thus Zuber's original expression gave values about 14% less than we 
would expect to obtain experimentally. 

The Problem of the Finite Flat Plate 


Equation (15) should give good results for finite plates as long as 
their dimensions greatly exceed \ d . However, as the plate dimen- 
sions approach ^ d[ _, A j/ A h will begin to vary discontinuously with plate 
size. The reason for this is that a small plate will have, say, one 
single jet on it. As the plate is increased in size, an additional jet 
(or jets) will be added, and the peak heat flux (which has been dropping 
up to this transition point) will suddenly be increased by a factor equal 
to the new number of jets. 

In counting the maximum number of jets that can be accomodated 
on a plate of a given size, we assume that, if there is more than one 
jet on a plate, then no jet will lie closer than Xq/4 to the vertical side- 
wall which must rise from the edge of the plate. Thus (with reference 
to Fig. 3) a square plate may be any size up to AX 2 and it will have only one 


jet on it. At this point q max will be very small — only a fourth of the 
predicted value. But now the single jet will give way to an array of 
four jets and will once again assume the predicted value . 




Fig. 3 Jet configurations on square heaters of various sizes. 


This rational can be continued. However, it becomes less satisfactory 
as larger and larger numbers of jets are considered, because the wavelength 
can vary somewhat around * d . The actual heat flux will equal the infinite 
flat plate value multiplied by (Nj \^/A h ), where Nj is the number of jets 
on the plate. This factor will approach unity as Nj and A h become large. 

The same rationale can easily be applied to circular, or other shaped, 
plates. 

Experimental Determination of q 

rn d/\ _ 


An important point relative to Zuber's equation is that it was never 
systematically tested against data obtained in the configuration for which it 
was intended. To approximate an infinite flat plate experimentally, one 
must first employ a very clean finite plate, much larger in any dimension 
than Xy. Then he must employ vertical side walls to prevent a horizontal 
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inflow of liquid since this has been shown [41, 4] to seriously influence 
Sri ax* 

The data we have located which meet these criteria are few. The 
vast majority of available flat plate data were obtained with strip or disc 
heaters in open pools and are hence unusable. The classical data of 
Cichelli and Bonilla [42] are for the correct configuration — a 9-1/2 cm 
dia. disc heater which formed the bottom of a tubular container for the 
boiled liquid. But very many of their data must be eliminated because 
they were obtained on "dirty" heaters. Most of the remainder are for 
nominal fluids of extremely low ourity — actually mixtures for which pro- 
perties are not known and correlations cannot be applied. Only a few of 
their data for ethanol remain for use. Berenson presented similar data 
for CC1 4 and n-pentane on 5 cm dia. heaters that were subject to very 
close control of surface condition. 

Owing to the paucity of available data, a program for the experimental 
determination of the peak heat fluxes on finite and "infinite" flat plates 
was carried out, both at earth-normal gravity and at elevated gravities. 

By "infinite" flat plates, we mean ones with dimensions much larger than' 
the "most susceptible" wavelength, and one with side walls to avoid 
any inflow of liquid from the sides. The flat plate we used accommodated 
more than three wavelengths for the alcohols, and about two for water, 
at earth-normal gravity; while these numbers increased as g 1 /^ at higher 
gravities. The observations were made for both viscous and nearly- 
inviscid liquids. The reagent grade liquids used were: acetone (CHg- 

COCHg), benzene (CgHg), isopropanol (CH 3 CHOH CH 3 ), methanol (CH 3 OH), 
distilled water and cyclohexanol (CH 2 (CH 2 ) 4 CHOH). 

A detailed description of the flat plate apparatus is given in reference 
[27], and we shall only describe some of its major features here. The 
apparatus consisted of a heater, a heater support, an emergency shut-off 
mechanism and a reflux condenser. This entire assembly was designed 
to operate either at earth-normal gravity or to be placed in the Centrifuge 
Facility of the Boiling and Phase-Change Laboratory at the University of 
Kentucky. A schematic diagram of the apparatus is shown in Fig. 4. 

The heater itself was made of pure copper with a circular bell 
shaped configuration. It has a diameter of 6.35 cm at the boiling surface, 
10.16 cm at the butt end, and a height of 10.8 cm. A pyrex glass cylin- 
der allowing a liquid head of about 6 cm sits on the heater surface and is 
held in place by two stainless steel flanges. Stainless steel has a much 
lower thermal conductivity than copper, hence the flanges have a very 
small cooling effect on the edges of the boiling surface. The total height 
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cock to control vacuum 
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of the heater with the glass cylinder is about 19 cm and it has a final 
weight of 5 kg. 

Figure 5 shows a sectioned view of the heater. Its lower portion 
consists of three rings wound with 4 meters of 1.62 mm dia. nichrome 
wire, capable of dissipating 6 kw of power. Five iron-constantan ther- 
mocouples are positioned as shown in the figure to monitor the tempera- 
ture in various parts of the heater body. The top two thermocouples 
were used to determine the heat flux to the boiling surface. A third 
thermocouple was used for a counter-check on the heat flux calculations 
and a fourth thermocouple to show whether steady state conditions are 
reached. 

The thermocouple in the base monitored the temperature for an 
electronic emergency shut-off system. This system cut off the power 
supply after the peak heat flux was reached. This was necessary because 
this transition effectively insulated the surface and permitted the tempera- 
ture in the rings to run away. This system is fairly complex and full 
details are given in reference [27], 

Elevated gravity tests were performed on the existing centrifuge 
facility. Complete details of the centrifuge design are given in reference 
[1]. Figure 6a shows the centrifuge designed and built to obtain 100 
times earth normal gravity at a speed of 360 rpm. When the flat plate 
heater is put on the centrifuge, its boiling surface is 77.7 cms from the 
center of rotation. 

Electrical connections to the centrifuge are made through slip rings 
attached to the drive shaft. Those include twelve circuits for power 
input, voltage measurement, and for thermocouples. The drive shaft also 
carries a tube which is connected to a vacuum pump through a stationary 
vacuum seal . 

A strobe light is placed on a bracket attached to the stand supporting 
the centrifuge. It is triggered by a photo-electric pick-off which senses 
light reflected from a shiny metal piece attached to the centrifuge arm. 
Thus each revolution of the arm triggers the light and the heater appears 
to stand still so one can view the boiling process. Figure 6b shows the 
flat plater heater assembly in position on the centrifuge arm. 

Before each test, the boiling surface was polished with 220 grit 
size emery paper to remove any rust or carbon deposit formed during 
the previous observation. The surface was then cleaned with soap and 
warm water and rinsed with the test liquid. This was done to ensure 
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thermocouples numbered 1-5 
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a.) Centrifuge facility — general view 



b.) 


Top view of flat-plate heater assembly 
installed at end of centrifuge arm. 


Fig. 6 Photographs of Apparatus 
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that there was no uncleanliness on the surface which might affect the 
wettability and affect the qp-^x data. The very careful experiments of 
Berenson [35] have shown that, in this configuration, is insensitive 

to surface conditions and that these precautions are more than ample. 

The resistance heating wire was energized and the current was 
gradually increased in steps. At each step the current and voltage were 
noted and thermocouple observations were recorded and updated until 
steady state was reached. It usually took about 5 minutes to reach 
steady state after the current was increased to the next higher value. 

Various regimes of nucleate boiling (e . g . the isolated bubble regime 
and the transition from single bubbles to slugs and columns) were easily 
discernable. The transition from nucleate to transitional boiling was 
identified by noting a slowdown of the boiling process and a sudden con- 
tinuous increase in the thermocouple reading. The thermocouple reading 
kept on increasing for a while, even after the power was removed. In 
most cases, the power was cut off after visual observation of the boiling 
transition, rather than by triggering of the automatic shut-off mechanism. 

The peak heat flux was computed from the observations of the top 
two thermocouples. Observation of the third thermocouple helped in 
making a counter check on the heat flux calculations. The maximum 
probable error of the peak heat flux is calculated in reference [27], It 
is found to be 7.7%. 

The peak heat flux observations were made for both viscous and 
inviscid liquids, both in the centrifige and at earth— normal gravity. A 
few of the observations were repeated to check for reproducibility of 
results. Each time a new observation was made the above-mentioned 
procedure for cleaning the boiling surface was used and fresh test liquid 
was employed. 

Discussion of Experimental Results 

The available qq-|axp data from both this and other investigations are 

plotted against L/A^ (or L'/2 tt/ 3 ) in Fig. 7. This abscissa scale gives 
the heater size in terms of the number of wavelengths that it will accom- 
odate. Data for circular heaters (L = diameter) and one point for a square 
heater (reference [41], L = width) are included in the figure. The origi- 
nal data and details of their reduction to dimensionless form are tabulated 
in Table I. 

Only those data for which the liquid was fairly inviscid are included 
in the figure. The question of viscous effects on q will be addressed 
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vertical side walls. 




Table I. Peak Heat Flux Data on Circular Flat Plate Heaters 


Liquid 

g /g e 

p 

(k Pa) 

M 

1 

i 

L’ 

-5 

q xIO 

max 

f W •' 

Smax 

2tt/cT 

^max^ 

Acetone 

1 

98.58 

599 

3.71 

3.94 

1 .18 


1 

98.93 

599 

3.71 

4. 10 

1 .22 


4.97 

23.99 

548 

7.69 

3.44 

1.14 


12.30 

25.37 

548 

12.07 

4. 19 

1.07 

Benzene 

4.97 

18.06 

384 

7.36 

2.98 

1 .21 


8.72 

24.88 

410 

9.85 

3.31 

1 .05 


17.5 

27.09 

413 

14.07 

4. 19 

1.07 

Iso- 

1 

22.75 

138 

3.68 

2.08 

0.79 

propanol T 


46.53 

182 

3.78 

3.34 

1 .0 



73.21 

233 

3.85 

3.82 

0.98 



88.52 

243 

3.88 

4.54 

1.10 



98.72 

254 

3.91 

4.26 

1 .0 


4.97 

18.06 

83 

8.2 

5.39 

1 .51 

Methanol 

1 

44.60 

370 

3.56 

3.69 

0.94 



47.43 

377 

3.56 

4.13 

1 .02 



88.04 

420 

3.72 

5.33 

1 .04 



96.38 

435 

3.75 

5.33 

1 .01 



96.52 

435 

3.75 

5.33 

1 .01 



98. 10 

439 

3.75 

5.33 

1 .0 

Distilled 

1 

14.48 

757 

2.24 

3.34 

0.70 

Water 


25.37 

930 

2.27 

4.45 

0.72 



29.51 

976 

2.28 

4.26 

0.66 



36.54 

1040 

2.29 

4.82 

0.67 



39. 16 

1060 

2.29 

4.54 

0.62 

• 


42.74 

1079 

2.30 

1 

-i 

■L - - 

4.54 

0.59 


While plotting the data for iso-propanol in Fig. 7, a small correction 
for the effect of liquid viscosity was made using Fig. 50. 
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in Chapter V, and some data for moderately viscous isopropanol are incl- 
uded but are corrected for viscosity in Fig. 7. Suffice it to say for 
the moment that a liquid is inviscid for this purpose when a viscosity 
parameter, M, is on the order of 400 or larger. 


With reference to Fig. 7, we see that all the data are very closely 
represented by our version of Zuber's flat-plate prediction, equation (15), 
when L Aq exceeds about 3. The two data for smaller LA d deviate 
sharply below the predicted value of q max = 1.14 q^^ > however. 


These deviant data — one is our own for a circular heater, the other 
was obtained years ago by Costello [41] for a square heater — these data 
serve strongly to support our present case. In Fig. 8 we give the peak 
heat flux for various numbers of jets using the expression suggested 
earlier. 


max 


N. X‘ 
J d 


(1.14 q, 


max- 


(16) 


Lines have been sketched in for various numbers of jets. No line has 
been carried back beyond the point at which the preceeding curve also 
rises above 1.14 q max . Likewise, no line has been carried forward 


beyond the point at which the following curve also falls below 1.14 q! 


'max 


We note that as LA d goes to about 3 and more, less variation in q max 
encountered . 


Z 

is 


The data point for LA d ~ 2.3, which actually represents 6 indivi- 

dual observations, corresponds very closely with one of the possible 
predictive curves. This success in rationalizing such a startling deviation 
from conventional expectation suggests that a more extensive comparison 
with experimental data would be warranted. 

Therefore, since we have only one other data point for small LA d , 
it is necessary to acquire additional experimental data. We turn to an 
unpublished study by D. M. Riherd [43], details of which we shall summa- 
rize briefly here. 

The apparatus is shown in Fig. S. It consists quite simply of a 
pair of electrodes positioned to accommodate nichrome resistance heater 
plates of varying width. Enough of each plate is short-circuited, as 
shown in the figure, so that the remaining active area is square. A 
chimney is fitted over the plate to contain the boiled liquid. Within the 
chimney are included spacers so the inside area of the chimney rises 
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heater width in wavelengths, L'/( 2 tt^/z ) 

Peak heat flux on a circular finite flat plate with 
vertical side wal Is. 









directly upward from the active heater area. 

Each plate was washed with soap and warm water after it was 
installed, and then rinsed with the test liquid. This was done before 
each set of observations was made on the plate. The results of these 
observations are presented in both raw and reduced form in Table II. 
The test liquids included reagent grade isopropanol and methanol, and 
95% pure ethanol. 

In Fig. 10 we have combined Riherd's data for square plates with 
Costello's point for a square plate. We also include equation (16) for 
Nj = 1,4,5,9,13, and 16. 

This figure shows very clearly the effect of finite heater size on 
q maX ' The data cling very closely to the line, between L/X^ = 1.0 and 
2. 0 .indicating that only one jet can be accommodated on the plate. This 
was substantiated by the scorch-marks on the used plates. Circular 
marks, i^/2 in diameter, were observed repeatedly. 

As the plate is shrunk to L = X^/2 we see that the data fall away 
from the curve . This is no surprise since the simple inviscid 
model no longer applies. The liquid return route is being squeezed into 
a viscous film crowded against the chimney, and considerations not 
included in the theory dictate q 

max 

This result provides a dramatic verification of the hydrodynamic 
theory. It reveals a 5-fold variation of q max in conformity with the 
theory. It is, perhaps, ironic that Costello used the data point shown 
in the figure to support inis own doubts as to the validity of the theory. 

Conclusions 


1.) The peak pool boiling heat flux on an infinite flat plate is 
given by 


max 


F 


1 14 t ... H I — — 

~ — : /c. h F _ mg (o P - o ) = I 

24 g f g v ^ f g ; 


14 


Sm 


ax. 


2.) If the plate is fewer than 3 wavelengths in width, the 
variations in Q rnax as jets are added or removed becomes very impor- 
tant. In this case: 

N . X ^2 

q max ~ (1.14 q maXf _) 


24 


Table II Riherd’s Observations for Finite Horizontal Square Plates 


Liquid 

width 

(cm) 

j Power 

* (watts) 

! JL_ 

! x d 1 

j- j 

q rnax 

Snax^ 

I 

Isopropanol 

CD 

5 

I 

156 

• 

CO 

fO 

0.777 



148 


0.738 


! 

156 

r j 

0.777 


1.52. 

102 

.93 i 

1.021 



105 


1 . 049 

Methanol 


132 

.89 

1 .048 



125 


.991 



127 


1 .010 



130 


1.029 

Ethanol 


130 

.90 

1 .207 



130 


1 .207 



130 


1.204 



119 


1.103 

Iso propanol 

1.84 

110 

1.12 " 

0.756 



108 


0.743 



111 


0.759 

Methanol 


146 

1.07 

0.789 



152 


0.825 



147 


0.800 

Isopropanol 

0.89 

74 

. 54 

2.183 



80 


2.354 



78 


2.292 

Methanol 


74 

.52 

1.726 



78 


1.817 



74 


1 .726 

Ethanol 


80 

. 33 

2.178 



82 


2.238 

Methanol 

1.21 

100 

.70 

1 . 266 



100 


1 .266 



100 


1.266 

Isopropanol 


101 

.73 

1 .614 



94 


1.506 



101 


1 .614 

Methanol 

2.16 

118 

1.26 

0.467 



1 14 


0.451 



110 


0.435 



116 


0.459 
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heater width in wavelengths L/X^=lI/( 2 ir/ 3 ) 

Fig. 10. Peak heat flux on a square finite flat plate 
with vertical side walls. 



III. THE PEAK AND MINIMUM POOL BOILING 
HEAT FLUXES ON FINITE HEATERS 


General Considerations 


The hydrodynamic method can be applied to the sort of finite heaters 
which are totally encompassed by a liquid bath, but some serious com- 
plications are introduced to the problem. References [5], [15], and 
[24] deal with such heaters and provide the basis of the present discus- 
sion. 


The present considerations all begin with a general peak heat flux 
expression obtained by combining equations (3), (5), and (1), namely 


Sri ax 
q 

max^ 


24 f 2tt 



(16) 


When X^ = X^ anc * ^j^h = this reduces to equation (15), however 

for finite heaters we must reopen the question as to what X^ and A/A^ 
are. These questions caused no problems in dealing with either the 
infinite or the bounded flat plates. 

Figure 11 shows the finite configurations under consideration here. 
They include cylinders, spheres, and horizontal ribbons oriented verti- 
cally with and without insulation on one side. Each of these introduces 
different ramifications to the problem. Much of this complication can be 
clarified by reviewing the derivation of q on the horizontal cylinder. 

Mi QA 

The Horizontal Cylinder, An Illustrative Example 

The geometry of a finite body exerts an important influence both on 
the size and the arrangement of the vapor jets. Photographic observa- 
tions of boiling on cylinders indicate — at least on the larger cylinders — 
that the jets adjust approximately to the width of the cylinder (plus the 
thickness of the vapor blankets, 2^) as shown in Fig. 11. If the cylinder 
is small the jets are small and the assumed spacing can reasonably be 
assumed equal to X^. As the cylinder increases in size the spacing must 
eventually spread to beyond X^ to accommodate jets which now exceed 
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X /2 in diameter . Sun showed that the spacing was 
diameters or about 4(R+6) in this case. Thus 


about two jet 



A h 


c R + 5 ) 2 

2R X 

d 


small 

cyls. 


A 

A h 


R+6 

8R 


large 


i cyls. 


( 17 ) 


Of course, equations (17) are true only insofar as Rj ~ R+6. 

Furthermore the wavelength, X^, is the dominant disturbance in 
the interface between the jets on large cylinders and it is picked up by 
the jets. The minimum Rayleigh disturbance, 2ttRj, is longer than X^ 
and would normally become Helmholtz unstable at a lower vapor velocity, 
U|_j. However, photographic evidence confirmed that vapor jets on large 
cylinders were much too short to have collapsed by virtue of the Ray- 
leigh disturbance. This means that the shorter waves of length, X^, are 
already well developed at the outset, while the Rayleigh waves require 
some distance to develop. Accordingly, Sun used X|_| = 2ttRj ^ 2n(R+6) 
for the small cylinders and X^ = X^ for the large ones. Using these 
X H 's and equation (17) in equation ( 16), gives 


q 

max 


^max-. 


6 

tt2 / 3 


(R’+ A) 3/2 
R 


and 


small 
cyls . 


3 3/4 


R' + A ! 

R’ j 

: large 
cyls. 


(18) 


where A : = 6 v/ g(p f - Pg)/a, a dimensionless blanket thickness. The transi- 

tion between small and large cylinders occurs somewhere in the neighbor- 
hood of X d = 4(R+6) or R' ^ 2.5 depending on the magnitude of 6. 

The parameter R’ , which has been variously named "the Laplace 
number," "the Rayleigh number," and the square root of "the Bond 
number", characterizes the ratio of buoyant forces to capillary forces 
in a system. As R* becomes very large, the system should approach 


4 The most susceptible wavelength used here is that for a flat interface. 
We do not include the correction for circumferential curvature [33,6] 
because that effect becomes lost in the irregularities of the interface 
near q 

max 
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a state in which it is no longer subject to capillary forces. In this 
state we would expect to see no further influence of R' upon q^.^/ 
^max z * This is what was four| d to be the case in [6], As R' became 

large. Sun measured 6^ 0.233 R so equation (17) gave Aj/A h — 0.155 
and c l max / c l maXz approached a constant value of 0.894. 

For small cylinders Sun approximated the measured values of 6 
with a fairly complicated equation in the form A = A(R'). Substitution 
of this expression in equation (18) gave: 


q , 

max 

= 0.89 + 2.27 exp (-3.44/^) (19) 

max 

Z ! small 
' cyls. [8] 

which represented approximately 900 data from a large variety of sources 
with very good accuracy. 

The Evaluation of A... and A. /A, 

H J h 


The fact that is the correct value of X^ for use with large cylin- 
ders can be used with any other large heater configuration since the 
situation in the upper interface is similar. Accordingly we can enter 
equation (16) with X^ = X^ and obtain in general: 


q : 

max 

i 

q 

max-, 

^ | large 

heater 


24 

TT i/3 



( 20 ) 


Likewise we may enter equation (16) with = 2 ttr. for any small 
heater. Then, using the fact that Rj must equal J ^A h / n )(A ./A h ), we 
obtain 


max 


max- 


4,- 
24 i 

~ TT ' 


na 


small 
' heater 


A h g(p f -P) Va h 


^ 3/4 


( 21 ) 


The problem of making general predictions of q,^-^ on finite heaters 
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then reduces to the problem of estimating A^/A h for any given situation. 

As a first step to determining Aj/A h we shall assume that the 
speed of the vapor passing through the blanket equals that in the esca- 
ping jet. For the speeds to differ would require the existence, both 
of pressure differences within the vapor escape path, and of significant 
dissipative mixing processes in the jet. We do not believe it is reason- 
able to took for either, and therefore assume that 6 simply adjusts to 
give equal velocities in both passages. 

For the large cylinder, this assumption, combined with a simple 
continuity statement (velocity times cross-sectional area is constant) 
gives 

2[4(R+S) 6] = A/2 (22) 

and for any sphere, large or small, it gives 

2tt(r+ 6/2 s ;6 = Aj/2 (23) 

For the small cylinder such a balance is not feasible since the 
vapor must flow horizontally in a long annulus subject to pressure drops. 
But for the small sphere equation (23) will still be true. From this 
point two paths can be followed: 

The path followed in reference [6] was to assume a jet configuration 
in terms of 6, and then to complete the derivation using observed values 
of 6. References [6] and [15] give the needed measurements of 6 for 
both cylinders and spheres as scaled from photographs. These data, 
and two additional points scaled from photographs in other papers 
[44,45,46] are combined in Fig. 12. Approximate lines have been fit 
through the data in both the large and small R' ranges. The results are: 


A small cylinders = (/3. 72/R' - 1) R' (24) 5 

A R 

large cylinders = 0.244 R' (25) 

^ small spheres = 0.20 R' (26) 

^ large spheres = 0.134 R' (27) 


5 

Sun used a more complex fit to the data — one which fit well in the mid- 
range but was very nearly equal to equation (24) for R'<_ 1. We are pre- 
sently more interested in low R' behavior than in transitional behavior at 
higher R' . 

® This result is a little higher than Sun's and represents a slightly better fit. 
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Data for water, nitrogen, and four organic liquids 
cylinders: • [6], -^measured from [44], 

(§) measured from [4 5] 
spheres: □ [I5 ],b measured from [46] 

0.244 R'- 

iyw-') 


.15 .2 



.3 .4 .6 .8 I 2 3 4 

Dimensionless radius, R' 


8 10 12 


Fig. 12 Vapor blanket thickness measured on spheres and 
cylinders. 
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Using these equations in equations (22) and (23) leads to: 


A 

j 

A. ; 
n | large 

I cyls. 


0. 155 


( 28 ) 


A. : 

— | = 0.143 (29) 

h | large 
spheres 


(30) 


The other approach is to combine the description of the assumed 
configuration of the jets with equations (22) and (23), solve the result 
for 5, and then obtain q max /q max from equation (18). This is reasonably 

safe to do in the case of large cylinders since the physical model 
(Rj = R + 6) is consistent with Sun's photographs. The result, 6 = 

0.244R, corresponds precisely with the experimental value given in 
equation (25) and leads to 


A. 

J 


= 0.220 

small 

spheres 


%ax 

q rnax F 

large 
cyls . 


0.904 


(31) 


This is negligibly higher than Sun's result of 0.894, and it is a com- 
pletely theoretical expression. While this result is accurate, the minor 
errors in the assumed characteristics of the vapor escape configuration 
accumulate more than we would like, in other cases. Such errors are 
particularly troublesome for the small heater configurations. 

Using the first approach, we combine the first of equations (17) 
with equation (24) and obtain 

A. 

77 - j =0.171 (32) 

h small 
cyls. 
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It appears from these Aj/A h results, that for all finite heaters 
except small bluff bodies (such as spheres), we can anticipate that: 


A. 

— - . 0.155 ± 10% 


Accordingly equation (20) gives 


q l 

%ax 7 

large 
heate r 


0.9 


( 33 ) 


(34) 


If we note that A^ for long slender heaters is ^ times the perime- 
ter, P, then equations (21) and (33) give 


q 

max 

q max-7 

long 

slender 

heater 


1.4 

Vp* 


(35) 


where P' = P/g(P f - P g )/ a 

p 

For small bluff heaters, A h ~ L , where L is an appropriate cha- 
racteristic length. Accordingly equation (21) gives 


max 

max 7 

^ small 

bluff 

heater 


constant 

A? 


(36) 


where the constant can be determined by the correlation of a few available 
data. If Aj/A h can be predicted, then there is no advantage to equation 
(35); a direct prediction can be obtained from equation (21). 

With these generalizations in hand, we now return to the specific 
geometrical configurations shown in Fig. 11. 
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Comparison of the Present Modified Theory With Data for Cylinders 


The present formulation will differ very slightly from that of Sun. 
The difference lies only in the fact that we do not use exactly the 
same description of his 6 data as he did. Using equation (32) in 
equation (21) we obtain 


q max 

Snax^ 


small 

cylinders 


0.94 / A S r? 


(37) 


(while if we had substituted P* = 2ttr’ into the simplified expression, 
equation (35), we would have obtained 0.885/V^R 1 — only 6% below the 
more accurate result). 

In Fig. 13 we compare equations (37) and (31) with the data given 
in reference £6]. The predictions are upheld perfectly. However, two 
general warnings are made evident in the figure: 

1.) Neither prediction should be applied beyond the point of where 
they join. This junction occurs a little beyond R' = 1 in this case, and 
at R' between 1 and 6 in all the cases we shall consider here. Thus, 
the transition generally occurs when the Taylor wavelength is comparable 
to the characteristic length of the heater. We also note that the sharp 
break between the two predictions represents the actual transition fairly 
well even though the actual transition is smooth. This will also prove to 
be true in the other cases we shall treat. 


2.) The prediction fails to be valid when R' is less than about 0.15. 
The reason for this will be treated at length in Chapter VII. Briefly, 
when R’ << 1, inertia proves to be no longer important in the peak heat 
flux mechanism. The hydrodynamic instability processes break down and 
the present predictions no longer apply. For that reason we claim no 
validity for any of the predictions when R* <0.1. 


Comparison of Theory with Data for Spheres 


For large spheres we can substitute equation (29) into equation (20) 
to obtain: 


max | 

q max z j large 

1 snheres 


0.84 


(38) 
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Dimensionless Peak Heat Flux q max /q max 


N 


5.0 


2.0 


10 


0 5 


r- 



t — i t~tttt [ 1 — i — rym 

range of about 900 data — 
half from [6] and the remain- 
der from 9 additional sources 
as reported in [6] 

—0.94 /(R') l/4 



equation (37) 
small R' 


^/^yp^777777 ^ 



equation ( 3IJ_ 
large R' 


0.2 


x i i I i ml i Li I i i l.i 


.904 


0 


0,*> I .5 

Dimensionless Radius R‘ 


10 


Fig 13 The peai< hua* flux on horizontal cylinders 
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(which is within about 6% of the simplified approximate expression, 
equation (34). 

For small spheres we substitute equation (30) and = 4TTR 2 into 
equation (21) and obtain: 


q 

max 

^max.^ 


small 

spheres 


1 .734 
/R 5 ' 


which is a special case of equation (36). 


(39) 


Since there were insufficient data in the literature to test these 
expressions adequately, Ded had to do his own set of experiments. We 
have already alluded to this work in the context of equations (26) and 
(27), and we should now describe it more fully. 

Figure 14 shows the general layout of the experiment used by Ded 
to measure q^^ and 6 in both nitrogen and double-distilled water. The 
test heaters were 0.635, 1.27, and 2.54 cm dia. copper spheres with 
ion^constantan thermocouples mounted in the center. The water-quenched 
spheres were heated to 510°C in the oven shown, and the nitrogen- 
quenched spheres were stabilized at room temperature. In both cases the 
spheres were then quenched in the test beaker. 

For boiling in water, a 12.7 cm I.D. glass beaker was placed on a 
hot-plate which kept the water saturated. A 15.2 cm I.D. double -walled, 
clear-glass, vacuum -insulated Dewar was used to hold the liquid nitrogen. 
It was kept saturated by heat leakage from the atmosphere, and the 
double wall prevented condensation on the outside so that clear photographs 
could be made. 


Details of the sphere mounting are shown in Fig. 15. An additional 
thermocouple was placed near the surface of the largest sphere for use 
in some of the water runs. This facilitated an experimental indication of 
the extent of the failure of the lumped-capacity (or uniform sphere temp- 
erature ) assumption. The worst case (in water) yielded only a 2% 
temperature difference between the two thermocouples. 

The spheres were polished with 500 grit emery cloth and the surface 
was cleaned carefully before each run. Since runs made with a sphere 
cleaned in acetone, and runs made with one simply sanded and then 
washed in water, gave identical results, both procedures were used. The 
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Fig. 14 Schematic diagram of the experiment. 




thermocouple leads 



Fig. 15 Test spheres for water and liquid nitrogen. 



acetone rinse procedure was used in obtaining most of the data available 
in the literature. 


The basic output of the experimental runs were temperature-time 
curves. The heat flux was obtained from these curves by noting that 


q = -p, 


in R 3 

4TTR 2 F 


(T) 21 = -p 
dt 


h 


(T) JL H_L 


dT 

dt 


(40) 


This method is called the "transient calorimeter" technique. Its 
use involves making the tacit assumptions that the Biot Number, Bi = 
hR/k w , is small and that the process is quasi-static. Bergles and 
Thompson [47], quoting from earlier sources, suggest that Bi at q 

rn 9x 

should not be much more than 0.4 for accuracy, and they warn against 
the second possible source of error — namely that the boiling process 
might not be quasi-static during the quench. 

A subsequent paper by Veres and Florschuetz [48] provided evi- 
dence that the present data were obtained under conditions for which 
the cooling can legitimately be assumed quasistatic 7 . When Bi exceeded 
0.4 in the present tests, equation (40) was abandoned in favor of the 
following computational method: 


The transient heat conduction equation for a spherically symmetrical 
body was solved numerically beginning with a known uniform temperature 
distribution at t = 0, and using the boundary conditions of zero tempera- 
ture gradient and T equal to a measured function of time, at the center. 
A simple Euler's method was used and the increments were chosen so 
that ccA t/ Ap 2 _ q. 139. Rerunning the program with larger and smaller 
values of ccAt/A r 3 did not alter the results. The instantaneous tempera- 
ture distribution was then averaged over the volume and T in equation 
(40) was replaced with this average value. An error analysis [14] 
gave 6.5 percent maximum probable error for the q __ data we present 

here. The maximum probable errors for q and R' were only 1.2 

rn 3.x r— - 

and 2.3 percent, respectively. I_ 


Ded's q data are presented in Table III. The temperature 
rn wX 

difference observed at q = q max is also reported for the reader's 
interest. With these data we include six points which he reduced from 


7 The rationalization of references [47] and [48] is somewhat intricate. 
It is discussed in [15]. 
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Table III. Peak Heat Flux Data for Spheres 


R 

(cm) 


q xlo 5 (W/m 2 ) 

max 

4T 

max 

(°C) 

! 

R ' | 

q 

max 

q 1 

max_, . 




Liquid Nitrogen, 

1 

, 2 

q =156, 000 W/m 

max 

— L- 

1 

| 

1 

.27* 



1 .18 

8.1 

12.01 

0.75 

1 

.27 



1 .25 

7.0 


0.80 

1 

.27* 



1 .39 

8.9 


0.885 

1 

.27 



1 .15 

9.5 


0.73 

1 

.27 



1 .06 

10.5 

» 

0.675 


.635* 



.95 

11.1 

6.0 

0.606 


.635 



1 . 14 

11.7 


0.73 


.635* 



1 . 12 

9.4 


0.71 


.635 



1 .23 

10.0 


0.784 


.635* 



1 .33 

9.2 


0.848 


.635 



1 .38 

10.6 


0.883 


.318* 



1 .43 

11 .7 

3.0 

0.91 


.318 



1 .79 

10.5 


1.12 


. 1 98 H 



1 .80 

8.9 

1 .87 

1 .14 


. 1 98 H 



2.50 

9.7 

! 

1 .59 


.0794* > 

H 


2.43 

7.0 

0.75 

1 .55 


.0794 h 



3.04 

9.2 

> 

2.14 


.0397*. 

H, 

T 

4.06 

7.8 

0.375 

2.60 


.0198 H) 

,T 


4.98 

7.5 

0.1875 

3.19 


2 

Water, q = 1 , 108,000 W/m 

max 


1 .27* 

i 8.20 

■ 1 

38.9 

5.06 

.739 

1 .27 

i 9.10 

43.9 


.82 

1 .27* 

10.09 

47.2 

> 

.91 

1 .27 

8.76 

40.5 

1 

.79 

.635 

7.76 

45.5 

2.53 

.69 

.635* 

8.93 

41 .6 

i 

.80 

i .635 

11.10 

40.0 

i 1 

i 

; 1 .01 

* denotes 

runs for which motion pictures 

were made 


H denotes 

results from 

t i me-te m pe ratu re 

data of Hendricks [19] 

T denotes 

tungsten-carbide sphere (all other spheres were 

copper) 
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temperature-time plots made by Hendricks and Baumeister in connection 
with a recent study of film boiling from small spheres [46] of both 
copper and tungsten carbide. 

These results have been nondimensionalized in accordance with 
equation (2), using property values as assembled by Reich and Lienhard 
[13]. The resulting q rna){ /q max values are plotted against R' in 

Fig. 16 along with the q max of Merte and Clark [49], Veres and 
Florschuetz. and Bradfield [50]. 

The predictions for spheres, equations (38) and (39), are also 
included in the figure. The comparison between theory and data is 
once again very good, indeed. 

The visual observations also serve to bear out the analytical 
models. Figure 17 shows typical photographs of boiling from spheres 
of increasing R' . The first picture. Fig. 17a, was taken from a 
motion picture supplied by Hendricks. The remaining three were 
obtained from Hycam movies (100 pps at 1/5000 sec exposure time.) 

In Fig. 17a one-half of a Rayleigh wave is apparent in the esca- 
ping jet. This wavelength is comparable with 2TTR as we anticipated. 
Figures 17c and 17d provide a side-view of the anticipated four-jet 
pattern for large R' although the jets are very blunt — much too short to 
contain a Rayleigh wave. Figure 17b is for a case close to the transi- 
tional range of R' and the single -jet configuration has deteriorated into 
an ambiguous conglomerate between the two models. These pictures 
are, of course, only static representations of processes that were 
repeatedly clear to eye when the movies were screened. The measure- 
ments of 6 mentioned on page 31 were also made from these movies. 

Comparison of Theory with Data for Horizontal Ribbons Oriented 
Vertically 


Consider next the case of a thin horizontal ribbon heater with 
the broad side of height, H, oriented vertically, as shown at the 
bottom of Fig. 11. We shall also give brief attention to such a 
ribbon with one side insulated. When H' is small equation (35) will 
apply, and the dimensionless perimeter will equal 2H' . When one side 
is insulated, P' = H' . Thus we obtain 

q max ; _ 1.18 1 

^max^ | small ^H* 

j horiz. ribbons 
* vert, orient. 


42 



Fig. 16 q on spheres. 






R =.0795 cm R 
liquid nitrogen 


R = .635 cm 
liquid nitrogen 


R = 2 . 54 cm R 
liquid nitrogen 


Photographs of vapor jets near q max on spheres of in 
creasing diameters in liquid nitrogen and water 




and 


max 


max. 


! small 

! horiz. ribbon 
! vert, orient. 

I 1 side insul . 


1 .4 

Vh t 


(42) 


For large h'’s we shall expect equation (34) to apply in either case. 

Table IV includes original data for vertically oriented ribbon 
heaters in four liquids: acetone, benzene, methanol, and isopropanol. 

The ribbons were all of nichrome, 0.023 cm thick, about 10 cm in 
length, and they varied in height, H, from 0.104 cm to 0.478 cm. 

They were operated as electrical resistance heaters and were connected 
to the heavy power supply electrodes through brass ribbon attachments 
which served to prevent vapor hangup by providing a smooth transition 
section. The range of H* ( = H / g(P f -P g )/ 0 ) was greatly increased by 
observing q max both at elevated gravity in the University of Kentucky 
Gravity -Boiling Centrifuge Facility, as well as at earth normal gravity. 


Complete details of the experimental method and apparatus can be 
obtained from [5] and [6], since exactly the same equipment and pro- 
cedure were employed. The probable experimental error in q^ax was 
about ± 4 percent although intrinsic variability of the data was as ± 15 
percent which is typical for such results. All ribbons had a smooth 
cold rolled finish (as delivered). Before each test the ribbons were 
carefully washed in soap and hot water, and then rinsed in the test 
fluid. 


We can be sure that, even on these small ribbons, q,^.^ did not 
occur prematurely by virtue of low thermal capacity effects such as 
Houchin [59] observed, since he was only able to observe the pheno- 
menon in water. Even though he used much thinner ribbons than we, 
he never witnessed the early burnout in organic liquids. 

These data are presented in dimensionless form in Fig. 18. 
Twelve additional cyclohexane data for this situation which were giv- 
en in [52] are also included in the figure. Equations (34) and (41) 
once again characterize the data perfectly. 


The peak heat flux was also measured on two 0.023 cm thick 
horizontal nichrome ribbons, vertically oriented, but with one side 
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Table IV Peak Heat Flux on a Vertically Oriented Horizontal Ribbon 


fluid 

H (cm) 

g/g 

e 

1 

! c 'max x1 ° 

H' 

q max ) 

9m ax ^ j 

Acetone 

0. 104 

i 

4.46 

± .16 

7 

; 0.65 

1.34 j 


0.130 

i 

4.90 

± .19 

j 0.80 

1.47 | 


0.207 

i 

3.70 

± .09 

j 1.27 

1.12 


0.354 

i 

3.36 

± .09 

2.21 

1.01 


0.480 

i 

3.07 

± .09 

' 2.98 

0.93 


0.366 

4.01 

4.42 

± .03 

4.56 

0.94 



8.30 

5.08 

± .06 

6.56 

0.92 


I 

17.84 

7.10 

± .03 

9.62 

1.03 


jL 

32.32 

8.40 

± . 13 

12.95 

1 .06 


) 

49.49 

9.06 

± .09 

16.03 

1 .08 

Benzene 

0. 104 

1 

3.96 

± .06 

0.64 

1.31 


0.140 

1 

3.70 

± . 13 

0.85 

1.23 


0.217 

1 

3.02 

± .09 

1.32 

1.00 


0.366 

4.01 

4.24 

± .06 

2.65 

0.99 



7.98 

5.18 

± .08 

3.70 

1 .02 



18.32 

6.36 

± .09 

5.60 

1 .02 


V 

31 .68 

7.45 

± .06 

7.35 

1 .03 


i 

49.49 

8.55 

± . 19 

9.25 

1.06 

Methanol 

0.140 

1 

6.32 

± .13 

0.88 

1 . 19 


0.252 

1 

5.05 

± . 09 

1 .59 

0.95 


0.354 

1 

4.40 

± .13 

2.25 

0.83 


0.480 

1 

3.90 

± . 13 

3.03 

0.74 


0.366 

4.01 

5.72 

± . 19 

4.63 

0.75 



7.98 

6.95 

± .05 

6.54 

0.78 



18.32 

8.86 

± .25 

9.91 

0.80 


V 

32.32 

9.85 

± .25 

13. 16 

0.77 

Isopropanol 

.0.354 

1 

3. 17 

± .09 

2.34 

0.76 


0.480 

1 

3. 10 

* . 13 

3. 16 

0.74 


0.366 

4.1 1 

4.81 

± . 16 

4.90 

0.81 



8.3 

5.66 

± . 22 

6.96 

0.80 



18.56 

6.60 

± .06 

10.40 

0.76 


i 

32.32 

8.04 

± . 19 

13.74 

0.81 
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obtained under elevated 
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heavily insulated with Sauereisen cement. A 0.252 cm ribbon was 
observed in methanol and a 0.478 cm ribbon was observed in acetone. 
The results were: 


at H' = 1.59 , 1.19 < ^?2 ax <1.30 

9 max ^ 

at H' = 2.98 , 1.03 — — < i .07 

'-'max ^ 

respectively. These data are plotted in Fig. 19 along with some high 
gravity data given by Adams [53] for higher values of H' in the same 
geometry. 

In this case we find that equations (34) and (42) give an excellent 
representation of the data. 

Both equations (41) and (42) correspond with Aj/A h = 0.155. The 
transition from large to small H' occurs at 6 when the ribbon is insu- 
lated and at 2.6 when it is not. 

CONCLUSIONS 

1. The method of hydrodynamic prediction of the peak heat flux 
on finite heaters is discussed in detail and certain general guidelines 
are set up for making such predictions. The assumption that the vapor 
velocities in the vapor blanket and in the jets must match, greatly 
streamlines prior descriptions and simplifies these guidelines. The 
guidelines can be summarized as follows: 


a.) 

q /q = 

max max-. 

CM 1 

TT 

0 

TT V 

A h 

9CP f - p g ) 

b.) 

q /q = 

max max-. 

24 

A. 

for large 

n ^3 

A h 

c.) 

^ 2 TT R for 

H J 

small 

he ate rs 

d.) 

V, = ^ for large heaters 
m d 

including 


A. 


3/4 


^ J 


for small heaters 


h 
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Fifl. 19 d m0K on horizontal ribbon heaters oriented 
vertically with one side insulated. 
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I 


e.) For small bluff bodies, A 


h 


L 2 and 


max 


max- 


_ constant 

small bluff ^ L ' 
bodies 


where the constant must be determined experimentally 
because Aj/A^ is only known in one configuration (small 
sphe res). 

f.) Except for the case of small spheres, A./A h seems to be a 
constant, very nearly independent of configuration. For the 
known cases (within about 10 percent) 


g.) Therefore, 


A. 


h 


£* 0. 155 


max i 


Smax 


~ 0.9 


Z : large 
heaters 


h.) q /q 

max max- 


_ 1 .4 


4 /R 


for long slender heaters. 


2.) Some special cases of the abo\& generalizations include the 
following predictions: 


a.) 

q /q = 

max ^max z 

0.94/ 4 /r' for small cylinders 

b.) 

q / q = 

max max z 

0.904 for large cylinders 

c.) 

q /q = 

1 .754/ v/R’ for small spheres 

max max z 

d.) 

q /q = 

max max z 

0.84 for large spheres 
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e.) 


W/%ax z = 1 - 18 /Vh- 

oriented vertically 


for small horizontal ribbons. 


f.) q /q = 1.4/ 4 /H’ for small horizontal ribbons 

max max^ 

oriented vertically with one side insulated 


g.) q /q = 0.90 for large horizontal ribbons oriented 

max max^r 

vertically, one side insulated or not. 

3.) The existing hydrodynamic qp-,^ predictions for finite bodies 
are summarized in Fig. 20. The figure includes an indication of the 
number of data by which each has been verified. This number exceeds 
the number of points actually shown in Fig. 18 and 19, since more 
than one observation has been lumped in some of the points. The 
curves have all been terminated at L* =0.1 on the left side since 
hydrodynamic predictions are known to deteriorate for <T 0.1. 
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Fig. 20. Collected predictions of q max for four finite heater configurations 



IV. A VISCOUS THEORY OF TAYLOR STABILITY 
WITH APPLICATION TO FILM BOILING 


Introduction 


We now turn to the question of the influence of viscosity in the 
hydrodynamic theory. As we indicated in Chapter II, a Taylor unstable 
process is common to both the peak and the minimum heat fluxes. 

Thus our study of the influence of viscosity in the hydrodynamic theory 
should begin with a viscous theory of Taylor stability. 

In 1S50, G.I. Taylor [54] discussed the instability of the horizontal 
interface between two ideal incompressible fluids of infinite depth. He 
showed that irregularities at the interface tended to grow if the accele- 
ration was directed from the less dense to the more dense medium. 
Bellman and Pennington' s contribution [31] was to extend his problem 
and showed how to take into account the interfacial surface tension 
and fluid viscosity. They gave closed form expressions only for the 
"most susceptible" wavelength (be., the wavelength for which the 
growth rate for a small disturbance at the interface is maximum) and 
the corresponding frequency. They obtained solutions for the situation 
in which only the surface tension at the interface was considered, but 
they were unable to solve the more general problem in which both the 
liquid viscosity and interfacial surface tension were considered. 

We shall begin with formulation of the instability problem similar 
to that made by Bellman and Pennington; we shall then incorporate the 
pressure contribution of curvature of a cylindrical surface as was 
done by Lienhard and Wong [33] for the inviscid case. Thus we shall 
attempt to obtain numerical values of the "most susceptible" frequency 
and corresponding wavelength as a function o' liquid viscosity, as well as 
radius of curvature of a cylindrical heater. Finally, we shall discuss 
the use of this evaluation in predicting the minimum heat flux. 

Hydrodynamic Analysis 

Before dwelling on the instability theory it is appropriate to state 
the various assumptions we are going to make: 

1) The fluids are incompressible 

2) The fluids are Newtonian 
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3) The fluid depths are large as compared to wavelength. This 
is not exactly valid for the vapor blanket thickness in film 
boiling, but we shall show later that the above assumption is 
not unrealistic. 

4) The upper fluid layer is a liquid and the lower fluid layer is 

a gas. We thus envision that the acceleration at the interface 
is always directed from the lighter to the heavier fluid. We 
assume this is so, because it is the situation encountered in 
all film boiling cases we will treat. 

5) Nonlinear effects are negligible. 

Figure 21 shows the interface between two incompressible viscous 
fluids of infinite depth. The linearized equations governing the motion 
are: 


u -} v — 0 
x y 


U t = 


1 u , 

p + — (u + u ) 

p x p XX yy 


1 . U, , . 

v = p ~ g h — (v +v ) 

t p y p XX yy 


These equations are satisfied by 


u = — < p - | 

x y 


and 


provided that 


and 



+ 

x 


p = p Q - pgy + p<p> t 


<p + <p =o 
xx yy 


-(♦ + t ) = I)! 

p ^ ' xx yy' T 


(43) 

(44) 

(45) 

(46) 

(47) 

(48) 

(49) 

(50) 
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r>c3 


Fig. 21 a. Interface between two incompressible 
viscous fluids of infinite depth. 


Fig. 2 1 b. Schematic diagram of a typica 
configuration of film boiling on 
cylinders. 



where u and v are the velocities in the x and y directions, respectively; 
p is the static pressure; p G is the total pressure, and p is the density 

of the fluid. 

Let us take, for the liquid phase: 



(p = A_ e ky+ lt cos kx 
f f 

(51) 


i' f = Bf e -m f y*'* 3 *- S in kx 

(52) 


p f = p Q - p f gy + 

(53) 

where k is the 

< 8 

wave number, & is the growth rate , and: 



m 2 = k 2 + P 'VP. 
f f f 

(54) 

Similarly 

for the gaseous phase: 



O = A e ky+xt cos kx 

g g 

(55) 


m n y+^t . 

ih = b e 3 sin kx 

g g 

(56) 

and 

P g = P o - p g 9y + 9 g<-Vt 

(57) 


m 2 = k 2 + P v/V- 

g g g 

(58) 


The real part of rrif or m g has to be positive so the velocity remains 
finite as y ± 00 . 


8 In Chapter II, we used to designate a true cyclic frequency. Thus 
the real part of exp (- Vt) was the cosine of the phase angle, and when 
exponential growth occured, was a pure imaginary number. Here we 
designate as '•*>, what we formerly wrote as (-i^). Henceforth w will 
designate a real growth rate — not a real frequency . 
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Considering that waves of height y = T]( X ,t) are propagated, we 
obtain for the linearized kinematic condition at the interface: 


•n 


t 



(59) 


or 


where 


d)t 

*n = v\ e cos kx 
o 


r, = + B f) 

O 


The boundary conditions at the interface, i.e. , at y=r|, are: 


u„ = u 

f g 


V = V 

f g 


-p 4- 2p (v ) = -p + 2 \JL (v ) - at) 

' r f y g 9gy xx 


(60) 

(61) 

(62) 

(63) 

(64) 


PfC(v f ) x -K Uf ) y ) = M g C(v g ) x + (u g y (65 ) 

The substitution of equations (51) through (53); (55) through (57) 
and (60) in equations (62) to (65) gives four linear and homogeneous 
equations in the four constants, A^, B P , A , and B • 

r f g g' 


pQCPf-Pgik 

L 'ju 



kA_ + m B_ ~ kA + m B = 0 

f f f g g g 


A + B_ + A - B = 0 

f f g g 

3 

2-, r g(P f -o )k Oh 

-o uu-2p k ;A + 2p km 

f f J f 0) m f 


+ [ p uu + 2p k 2 l A - [2 p km jB = 0 

g g g g g g 



( 66 ) 

(67) 


( 68 ) 
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and 




<- c O 
3 O P 
° T3 i 
n .E -O 


.S c 
^ <u 
0 ) 5 > 

C c w 
® O c 
£ n 5 s 


C ^ a) 

03 ■**- r 

5 o c 

C O) 

r c <d 


“0^0) 

0) (/) c 

N 'H ^ 

> to 5 

c (d ^ 
CO C r 

8 « c 

u - i I o 

X 4-1 

.21 • TJ 

St.® 

X „<D * 


i § * 


As proposed earlier by Lienhard and Wong [ 33 ], a two dimensional 
model can be used for cylindrical geometries. The effect of transverse 
curvature is treated in the form of an additional oscillating pressure 
difference component across the interface, as shown in Figure 22 . 

The expression for the transverse pressure may easily be written 



2R 2 ( 1 +d g /R) 2 


where d g is the depth of the vapor blanket over the wire 


(71) 



P 

g 


Effect of transverse curvature 
considered as a sinusoidally 
varying pressure. 

Figure 22 Contribution of transverse curvature on a two 
dimensional wave 
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The evaluation of the determinant equation, (70a) gives: 


r 3 

j_-g(p f -p )k + ok 


~~2 +(p f +P g )Uj2 J [^ f ( k + m f )+ V k + m g ) ] 


+ 4wk0a k +[I m l[[j k + iim] = O (70b) 

r 9 9 9 f f 

Our next step is to put equation (70b) in a more usable form, 
so that it can be solved explicitly for frequency and wave number, 
or wavelength. The vapor viscosity M g is much less than P f so it 
may be neglected in comparison to liquid viscosity. Thus we may 
write, for the growth rate, w , 


tU " ffp r p g’ P f — p g * M F>° • k > 9 , R c ) ( 72 ) 

This expression relates eight quantities which are expressible in 
three dimensions. Using the Buckingham Pi -Theorem, we can recast 
this problem in terms of five dimensionless groups. For these 
groups we choose: 


(i) a dimensionless growth rate or "frequency" 


0 = 0) ! 


r a 


.1/4 


L 3 

g (P f -P g ) 


(ii) a dimensionless wave number 

. 1/2 

'f V 


K « k r^-rl 

'.g(p -p J j 


(73) 


(74a) 


or a dimensionless wavelength. A, defined as 


1/(/3 K) 


(74b) 


(iii) a dimensionless liquid viscosity parameter 


M 


o a 


3/4 


V / 4 cp f - y 3/4 


( 75 ) 


The square of this group is very nearly the 
Borishanski number, N, as defined in equation (11). 


(iv) a dimensionless density 


Pf- p g 


(v) a non-dimensional cylinder radius 


, _ - g(p f - p 0 ) i 

R * = R , *- 

C C L. 0 J 


1/2 


(76) 


(77) 


This number is related to the Bond number. 
Bo, by Bo = R ' c 2 * 

Using these dimensionless numbers, we may write 
equation (70b) as: 


1 - K 2 + rA — -Att + 


K 


K 


2 Bo T K (K 2 +qm) 1/2 (K 2 +T,M) 1y/2 


Q 


40 K 


+ 


K 


2 1/2 2 1/2 
T(K + QM) 7 MF 2 Bo(K + 0M) ' 


= 0 


(78) 
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Equation (78) gives a relation between growth rate and wave 
number for an unstable disturbance (corresponding to il positive and 
real), when the interfacial surface tension, transverse curvature and 
liquid viscosity are taken into account. 

Prom equation (78) it is clear that when the wave propagation 
velocity is zero, there is no effect of liquid viscosity on the critical 
wavelength. Also, when M ”* 00 ( i .e . , the liquid is inviscid) equation 
(78) reduces to the same equation as obtained by Lienhard and Wong 
[33] for the inviscid case. 

We are interested in the "most susceptible frequency", or the 

value of uj for which growth rate of the disturbance is maximum (i.e . , 

the one for which !iL_ = 0). Differentiating equation (78) with respect to 

dk ^ 

k and setting _ = 0, gives 
dk 


-2K + -2-r + 1 


K 


3K 


TK“ (K 2 +fiM) 1/2 (K 2 +QM) 3/2 (K 2 +ftM) 1/2 


K 


^ + 


q 2 k 


4Q 1 


(K + fiM) / r(K 2 +QM) 3//2 Mr 2 Bo(K 2 +QM) 1//2 


K 


2Bo(K 2 +QM) 3/2 


= 0 


(79) 


Thus we have two nonlinear equations (78) and (79) in two unknowns, 
and K. The two equations were solved numerically on an IBM/360 
computer using a subroutine: XFNLES, which solves a system of 

nonlinear equations. This subroutine is available in the Numerical 
Analysis Library [57] for the S/360. 

For the physical cases that we shall consider subsequently, T is 
approximately 0.9995. This is the value that we used in the calcula- 
tions. We might just as well have set F= 1 , however, since r would 
have to be much farther from unity to alter the computations noticably. 
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The "most susceptible frequency" and corresponding wavelength 
are plotted in Figs. 23 and 24, respectively, as a function of the 
viscosity parameter M and Bond number. Bo. As evident from 
these figures, the effect of liquid viscosity is to increase the wave- 
length and to decrease the corresponding frequency. 

Figures 25 and 26 show the effect of Bond number and liquid 
viscosity separately on the frequency-wavelength relationship equation 
(78). Transverse curvature of the cylindrical heater reduces the 
"most susceptible" wavelength and increases the frequency. Figure 
26 shows that apart from increasing the wavelength, viscosity also 
tends to increase the region of near neutral stability slightly. By 
"region of near neutral stability" we mean the range of wavelengths 
that can exist within any specified range of frequency close to the 
maximum frequency. 

Experimental Determination of Vapor Blanket Thickness, Wavelength, 
and Growth Rate 

An experimental program was carried out to observe the wave- 
length, its rate of growth, and the thickness of the vapor blanket 
surrounding the wire heater during film boiling in viscous liquids. 
Since the viscosity of most of the liquids is fairly low when they 
boil at normal pressures, the experiments had to be performed at 
very low pressures to illuminate significant viscous effects. 

With the limitations of attaining maximum vacuum in the existing 
boiling apparatus in mind, a survey of various available chemicals 
was made. It was found that Cyclohexanol , CH 2 (CH 2)4 CHOOH, and 
Propylene Glycol, CHg(CHOH) CHgOH, were best suited for the 
purpose. However, nearly all the experiments reported here were 
performed with cyclohexanol as the boiling liquid. Complete infor- 
mation regarding the relevant physical properties of cyclohexanol are 
presented in reference [27], 

The basic element of the apparatus was a test capsule 8.9 cm 
wide, 8.9 cm high and 17.8 cm long, with glass windows in the 8.9 
by 17.8 cm walls. It was made of brass and insulated on the sides 
with 0.63 cm thick styrofoam sheet. An electric preheater and 0.63 
cm diameter brass holders to support the test heaters were fitted to 
the capsule. Nichrome wires were used as test heaters and copper 
hooks were attached to their ends so they could be mounted in the 
brass holders. The hooks were attached in such a manner as to 
minimize the end effects. A 2.54 cm marker was mounted on the 
bottom of the capsule to provide a reference dimension for the 
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reduction of photographic data. 


Figure 27 shows a schematic diagram of the apparatus. AC 
power was employed in most of the experiments. The power supply 
to the wire was calculated from the measured current and voltage in 
the wire. A mercury manometer or a vacuum gage was used to note 
pressure inside the capsule. This pressure was corrected to take 
into account the head of liquid above the wire. An identical apparatus 
was used by Lienhard and Sun to make similar measurements. Full 
details are given in reference [5]. 

Nichrome wires, about 10cm long, were cleaned with soap and 
hot water to remove any grease or oily matter and then rinsed with 
the test liquid. The wire surfaces were smooth and had a cold-rolled 
finish (as supplied by the manufacturer). The wires were mounted 
in the brass holders and the glass windows were positioned. The cap- 
sule was filled with reagent grade test liquid to a level of about 2.5 cm 
above the wire. The vacuum pump was started and the preheater was 
turned on to heat the liquid to saturation temperature. The preheater 
was turned off before energizing the wire. This was done to avoid any 
effects of convective currents and electric fields. 

The current in the wire was steadily increased until the peak heat 
flux was reached and the transition from nucleate to film boiling was 
observed. Thereafter the current was reduced until film boiling 
started to disappear at the ends. This insured that the heat flux was 
close to minimum. 

Still pictures of the film boiling phenomenon were taken and 
observations of the temperature and pressure were also made. The 
heat flux was again increased in steps. At each step a photograph was 
taken and other corresponding observations were made. This procedure 
was continued until the heat flux was such that the wire started glowing 
or the bubble merger mechanism was evident. By "bubble merger" 
we mean that two neighboring bubbles growing in the first half of the 
cycle come into contact with the adjoining bubble from the second half 
of the cycle. Again the heat flux was decreased in steps and the 
process repeated. 

This procedure was repeated for wires of different sizes and 
various pressures. Each time a new wire was used, the liquid in 
the capsule was also replaced. In some cases high speed movies were 
made, to facilitate study of the growth rate of the disturbance, with a 
Hycam motion picture camera at a speed of 500 frames per second and 
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Fig. 27 Schematic diagram of the apparatus. 




a n exposure time of 1/1250 second. 


Photographic information was used to make three different kinds 
of observations. We describe them under three separate subtitles: 

Wavelength measurements — Still pictures were taken with a 
2.54 cm marker as a reference dimension, so the wavelengths 
could be scaled from them . Care was taken to avoid situations 
where bubble merger or deformation of the wave pattern had taken 
place. This deformation may be caused by drift in the phase angle 
along the wire or by oscillations generated through the movement of the 
outgoing bubbles. The probable error in the measurements of the 
shortest wavelengths was only about ±4.5 percent. The error was 
less for larger wavelengths. 

Vapor blanket thickness measurements — In film boiling, a 
vapor blanket of finite thickness always surrounds the heater and there 
is no liquid contact with the surface of the heater. The instability of 
the upper boundary of this vapor blanket causes a steady release of 
bubbles. For the application of the theory developed in the previous 
section we must obtain the corrected radius of the heater wire by adding 
the minimum blanket thickness to R. 

To observe the vapor blanket thickness, representative pictures of 
film boiling were enlarged. The minimum diameter of the vapor blanket 
surrounding the wire was measured, using the 2.54 cm marker again as 
a reference. This was done to avoid any refraction effects that might 
creep in if the wire diameter were chosen as a reference. Furthermore, 
the wire diameter could not be seen clearly through the vapor blanket. 

The wire diameter was subtracted from this measurement to give 
twice the vapor blanket thickness. In other words, we anticipated that 
vapor surrounded the wire equally on the top and the bottom. Some- 
times a cusp appeared on the blanket in the wake of a departing bubble, 
and we have ignored it. 

Measurements of growth rate of disturbance — Hycam movies were 
viewed on a microfilm viewer. At first, a preliminary survey was made 
to sort out those bubbles which grew undisturbed by other adjoining 
bubbles. Starting with a frame in which a bubble had just broken away 
from the interface, the height of the interface was measured from the 
lowest boundary of the vapor blanket. Later, the minimum diameter of 
the vapor blanket tube surrounding the wire was subtracted from each 
of the above observations to give the amplitude of the wave. 
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The minimum height of the intenfa.ce was used as a reference 
dimension for obtaining dimensionless amplitude. Thus, at a parti- 
cular value of x, say x=0, one may write 


■n = T1 exp (UO t) 


or 


T| _ it. 4 

*\r 


= exp [Ot / g3 ( Pf - p )/ct] 


4 

but we can easily show that /p^ Pf - p ya = 1.612 ^ d , 


so 


(80) 


(81) 


In (JL )= 0(1.612 w d t) ( 81a) 

% F 

where r] Q is the minimum depth, d , of the vapor and «J d is the "most 

susceptible" frequency (or growthrate) for the disturbance in the invis- 
cid fluid in the absence of any geometrical effects of the heater. In all 
cases the dimensionless amplitude, Tl/r^, was plotted against dimen- 
sionless time, on semi-logarithmic graph paper. The slope of 

the curve at any instant gave the dimensionless frequency 0. The pro- 
bable error in the observation of the linear growth rate was ± 10 percent. 

Now we would like to compare the experimental observations of 
wavelength and frequency with the theoretical predictions. However be- 
fore we can do this it is necessary to present some sort of correlation 
for vapor blanket thickness, because the wire radius has to be corrected 
for it. 

Vapor Blanket Thickness Correlation 

Baumeister and Hamill [58], while analyzing heat transfer from 
wires in film boiling, developed an expression for vapor blanket thick- 
ness. Their theoretical model for the film boiling configuration was a 
fairly approximate one; a sequence of spherical domes connected by an- 
nular passages. Neglecting inertia, they solved the equations of motion 
and energy with the assumption that the heat transfer rate is maximum. 

For the heat transfer coefficient for film boiling on wires they 

got: 
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(82) 
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where A H dg/R. All the physical properties of the vapor in the above 
equation are evaluated at the mean film temperature, i ,e . , at T= T sat 
+ At/ 2. The latent heat of vaporization has been corrected to take 
into account any sensible heat required to heat the vapor above its 
saturation temperature. Assuming a parabolic velocity profile and 
mean film temperature for vapor in the vapor blanket, Breen and 
Westwater [59] gave a corrected expression for hy : 


h* = h 

fg fg 


i + 


0.34 c v A T -.2 


fg 


Baumeister and Hamill changed the constant in equation (82) from 
0.23 to 0.35 for small wires (R 1 < 0.05) to take into account any heat 
transfer from the domes. This constant was further corrected to 
0.485 for larger wires (R' s 0.05) to fit the experimental data. Essen- 
tially, the constant in equation (61) was fixed from experimental results. 
This was necessary owing to their necessarily simplified model of the 
film boiling process. 

Their theoretical expression for the vapor blanket thickness was 


A = exp 4 . 35j- 


k U AT 1/4 

g g 1 

h! 2R 0 o J 

f g g 


r 1 

i 2 

L 1 +9R' (1 +A) 


J/4. 


2 J 


}- 


(83) 


where the term in the second square bracket is close to unity for 
small values of R' and A. 

If heat flux, q, is taken to be an independent variable; the temp- 
erature difference. At, between the heater wall and the surrounding 
liquid can be written as q/h w so equations (82) and (83) yield: 
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The group ^ g /h*^ g p g a wea ^ function of temperature. We took 
A T= 425 °K, which is typical of the film boiling situation, and plotted 
the present data for cyclohexanol , as well as some old data of Sun 
for acetone [1] against equation (84). The predictions were too high 
and did not correlate the data well for various dimensionless radii. 

This suggests that we should take another look at the various factors 
influencing vapor blanket thickness. 

The vapor blanket thickness will primarily depend on seven inde- 
pendent variables, q, p g , h* fg , g(p f -p g ), a, p g and R. The heat flux 
q is imagined to be transferred by conduction and used fully in the phase 
transformation . The eight variables, including vapor blanket thickness, 
d g , can be written in four dimensions: joule, kg, meter, and sec. 

Thus, in accordance with Buckingham Pi Theorem, we expect the problem 
to be reducible to a relation among four dimensionless groups. For 
these we choose: 


a) The dimensionless vapor blanket thickness, A. 

b) The Bond number. Bo, or ratio of buoyant to surface tension 
forces . 


c) A ratio of viscous to surface tension forces. 


B 


q4 


g 
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(85) 


d) A ratio of inertial to viscous forces — a kind of Reynolds 
number based on the vapor velocity, 

Re s , ' q ^» (86) 

U h* v } 

g fg 

Baumeister and Hamill's equation (84) involved the first three 
groups but the fourth group was missing because they neglected interial 
terms in their equation of motion. Their assumption was realistic as 
long as either the wire radius, or the heat flux, or the vapor blanket 
thickness is small; but the ratio of inertial to viscous forces might be of 
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the order of magnitude of unity when the heat flux is high or the vapor 
blanket thickness is large. 

Radiation losses can also form a significant portion of the total 
heat transfer from a wire in film boiling. We employed nichrome wires 
in all the experiments performed. The total hemispherical emittance 
for slightly oxidized nichrome at a temperature of 750°K has been 
given by Touloukian [60] as 0.65. Thus the radiative heat flux may be 
written as: 


% - °- 65 a < T W 4 - T sa/> < 87 > 

For At ci 425°K, which is typical for our wires, equation (87) gives 
0.132 x 10 5 watts/meter^. This heat flux may form as much as 10-20 
percent of the total heat flux, and it is absorbed by the liquid. 

To know the depth of penetration of this radiant heat flux, absorp- 
tive properties of liquid cyclohexanol are needed. However, this infor- 
mation is not available in the existing literature. Values of the absorp- 
tion or extinction coefficient for water [61] indicate that all the radiant 
energy is absorbed in a very thin layer of liquid. Even if absorbtion 
is considerably less for cyclohexanol, the radiant energy absorbed by 
the liquid adjoining the interface will largely find its way back to the 
interface where it will be utilized in generating additional amounts of 
vapor. Thus we incorporate all of q r in q for use in equation (84). 

Depending upon the total amount of vapor generated, the vapor 
blanket thickness and vapor velocity in the vapor blanket, will adjust 
themselves in such a way as to optimize the process. The Reynolds 
number, Re, as defined in equation (86) incorporates the vapor velocity 
and the vapor blanket thickness, hence it is an important correlation 
parameter . 

In Figure 28, we have plotted 20 data points for vapor blanket 
thickness on wires in the range 0.11 < 0.65 during the film boiling 
of acetone and cyclohexanol. While the abscissa is the same as 
Baumeister and Harrnll's, the vapor blanket thickness has been corre- 
lated with the Reynolds number as an additional parameter. This 
correction has been obtained by trial and error procedure and as yet, 
no theoretical reason is offered for taking the exponent in the inertia 
correction term to be 2/3. In forming the vapor blanket thickness 
correlation, all of the vapor properties have been evaluated at the mean 
film temperature. The raw data are presented in Table V. 
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Table V Vapor Blanket Thickness Data for Horizontal 
Cylindrical Heaters at Earth Normal Gravity 


R 

(mm) 

R' 

qx 10 * 

x m 2 J 

-IV * 10 - 

0 h " O ' 

9 fcj 

_C]R 

f ; h ' 

9 fg 

Observed Vapor Blanket 
Thickness 

(minimum - maximum) 

ci ;-( rnrn ) 

i = d c /R 

Cyclohexane! 

0.2667 

0.144 

0.83 

10 

1.25 

0 . 12 - 0. 15 

0 . 45 - 0.55 

0.4128 

0.216 

1 . 16 

64 

3.42 

0 . 12 - 0. 16 

0 . 28 - 0.38 





1.70 

104 

5.55 

0 . 18 - 0.21 

0 . 44 - 0.51 





2 . 12 

130 

6 * 86 : 

0 . 19 - 0.23 

0 . 47 - 0.56 



0.227 j 

i .01 

1 1 .6 

2.30 

0 . 17 - 0.24 

0 . 40 - 0.58 

0.5144 

0.268 : 

O'. 96 

i 

53 

3.60 

0 . 13 - 0. 17 

0 . 26 - 0.34 




[ I 

1 .68 

92 

6.00 

0 . 21 - 0. 26 

0 . 41 - 0.51 



0.278 

0.67 

7.8 

1.90 

0 . 15 - 0. 19 

0 . 29 - 0. 37 



0.2 

87 

0.79 

3.4 

2.64 

0 . 17 - 0.22 

0 . 33 - 0.43 



0.294 

0.79 

1.8 

2.54 ! 

0 . 1 1 - 0.13 

0 . 22 - 0.26 



1 


1 .02 

2.4 

3.30 

0 . 10 - 0.13 

0 . 20 - 0.26 

0.6540 

0.343 

0.84 

46 

3.86 

0 . 13 - 0. 17 

0 . 20-0 . 26 

1 

l 


! 

0.92 

i 

50 

4 . 25 . 

0 . 21 - 0.26 

0 . 32 - 0.40 


J.. 

Acetone ' 


0.2667 

0.3200 

0.4128 

0.6440 

0.6540 

0.8250 

1.0280 


0 . 161 

1.04 

0.33 

1 . 44 

0.203 

1.30 

0.45 

2 . 27 

0.257 

1 .19 

0.42 

2.63 

0.409 

1.01 

0.35 

3.55 

0.458 

0.88 

0.31 

3 . 10 

0.514 

1.17 

0.41 

5 . 16 

0.650 

1.11 

j 

0.39 

6.20 


0 . 04 - 0.05 
0 . 04 - 0.03 
0 . 07 - 0.00 
0 . 05 — 0 . 06 
0 . 07 — 0. 10 
0 . 03 - 0. 1 1 
0 . 03 - 0.09 


0 . 15 - 0. 19 
0 . 13 - 0. 19 
0 . 17 - 0.18 
0 . 08 - 0. 1 
0 . 1 1 - 0. 16 
0 . 11 - 0. 14 
0 . 06 - 0.09 


+ 


Data reduced from observations of Sun as reported by Liennard 
and Carter. [1], 


The group, qRA/Hgh* f , varied approximately from 0.2 to 2.5 
in all the observations. Smaller values of this group usually occur 
near the left hand side of the figure. The data scatter represents 
minimum and maximum measurements of vapor blanket thickness. 

The accuracy of these measurements is ± 10 percent. The data seem 
to be correlated well by the solid line whose governing equation is: 
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( 88 ) 


The constant, 4.60, in equation (88) would correspond to a value of 
0.32 instead of 0.23 for the constant in front of equation (82). 

Comparison of Wavelength Predictions with Experiment 

Observed wavelengths for cyclohexanol boiling at a temperature of 
302. 5°K, corresponding to absolute pressure of 0.296 kPa, are displayed 
in Figure 29. The value of the liquid viscosity parameter, M, is 
close to 5. The theoretical prediction for the viscous and inviscid 
cases is also shown in the -figure. The wire radius has been corrected 
R c to take into account the thickness of the vapor blanket surrounding 
it. The raw wavelength data are presented in Table VI. 

Although wavelengths were measured for various heat fluxes, 
those displayed in Figure 29 correspond to the lowest heat flux. This 
is done to avoid the longer wavelengths which may, as we will see 
shortly, be favored at higher heat fluxes. The data show a wide 
variability, but the lowest points in the range of the data scatter do 
embrace the theoretical prediction which takes the dominant wavelength 
to be the "most susceptible" one. The inviscid predictions would have 
suggested wavelengths about 28 percent too short. 

Figure 30 shows photographs of film boiling for two of the data 
points in Figure 29. Although there is slight phase difference across 
the length of the wire, the Taylor wave is very well developed. The 
bubble release pattern is good. The process is slower and viscosity 
seems to have dampened the small interfacial disturbances which may 
be observed in less viscous liquids. 
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Fig. 29. Wavelengths on horizontal cylinders. 




Table VI Wavelength Data of Cyclohexanol for Horizontal 
Cylindrical Heaters at Earth Normal Gravity 
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Wavelength data for M = 16 are plotted in Figure 31. Here the 
absolute pressure is 1.06 kPa and the saturation temperature is 329°K 
In this case, observed wavelengths are higher and the data again show 
wide variability. This variability in data could be explained from the 
dispersion relation shown in Figure 25. Near the maximum frequency 
there is a wide region of near-neutral stability. Thus, for frequencies 
slightly less than the maximum, a large range of wavelengths close to 
^ is possible. 

Comparison of Wave Growth Rate Predictions with Experiment 

Figures 32, 33 and 34 display plots of dimensionless wave ampli- 
tude versus time on semi-logarithmic coordinates in an increasing 
order of R' . In all cases, the liquid viscosity parameter, M, for 
cyclohexanol is 16. The data contained in each figure is for two 
to three randomly-picked, regularly-growing waves during 5 to 10 
seconds of motion pictures of film boiling. 

All the figures reveal that bubbles grow linearly during the first 
12 percent or so of growth. This is the period during which we would 
expect our theoretical predictions of the frequency to be valid. Occa- 
sionally during the early growth of a bubble, an oscillation of the 
interface was observed. This was probably caused by the superposition 
of a disturbance left in the wake of a departing bubble. Such a 
disturbance is evident in the wake of the left-hand bubble in the bottom 
picture in Figure 30. Data corresponding to such disturbances have 
been identified with dotted symbols in the wave growth rate diagrams. 
One may also note that the bubble grows in height to about 75 or 80 
percent of the wavelength before it leaves the interface. 

Having measured the frequency corresponding to a linear growth 
rate and the dominant wavelength, we next wish to see how well the 
two compare with the predicted dispersion relation. In Figures 35, 36 
and 37 we trace the dispersion relations for the three dimensionless 
corrected radii and M=16, and we display the experimental points on 
them. The relation between wavelength and frequency is borne out 
qjite well in each case. It is clear that wavelengths with frequencies 
slightly less than the "most susceptible" frequency can easily occur. 

In Figure 38, observed frequencies are plotted as a function of R c ' 
along with the viscous and inviscid predictions. The data fall slightly 
below the viscous predictions, but the inviscid theory would have pre- 
dicted still higher frequencies. 
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Fig. 32 Growth of waves on a 0.406mm dia 

wire heater in cyclohexanol. P = l.06kPa, 
q = 0.85 x | O 5 watt/meter 2 , f b = 23 bubble/ 
sec., =45.5 hertz, R' = O.II, Rl=O.I85. 
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Fig. 53. Growth of waves on a 0.51mm dia wire 
heater in cyclohexanol. P= 1.06k Pd, 
q = 0.86 x|0 5 watt/meter 2 , f b =20 
bubble/sec, w d =45.5 hertz, R , = 0.I4,R' C =0.22 
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Fig. 34 Growth of waves on a 0.825 mm dia 

wire heater in cyclohexanol. P = 1 .06 k Pa, 
q = 1 .01 x 10 5 watt /meter 2 , f b = 18 bubble/ 
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cyclohexanol at 1.06 kPa. q=0.85x!0 5 watt /meter 2 . 
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Fig. 38. Experimental observation of frequency of dominant 
wavelength on cylinders. 




The inviscid theory shown tn Fig. 38 can be obtained by reducing 
equations (78) and (79) to the Lienhard-Wong relations. This in turn 
can be done by letting M " °°. Figs. 23, 24, and 26 show that when 
M > 50 viscous effects cease to be important. The results are 

1 - K 2 + 1/2 R' c 2 = Q 2 /r K (78a) 

and at the maximum fl, 

2K = fi 2 /TK 2 (79a) 

The desired relation between R' and for low viscosity fluids is 
obtained by eliminating K between these two equations. 

Figures 39 through 42 show the growth of waves in low viscosity 
fluids. Figures 39 and 40 for M = 52 and 68, respectively, were 
obtained from the present motion pictures while Figs. 41 and 42, 

(M = 264) were reduced from a movie by Lienhard and Wong [62]. The 
growth rates from these figures are compared with the inviscid theory 
in Fig. 43. The comparison is very satisfactory. That the growth 
rates do not exceed n max is encouraging, since we have already seen 
that 0 should be equal to or slightly less than its maximum value. 

As a matter of interest. Fig. 42 shows an oscillatory behavior of 
the interface during early growth, probably due to supe rim position of a 
disturbance caused by the preceding departing bubble. We have 
attempted to ignore the disturbance again and to draw the line through 
the average growth curve. 

From the previous observations of growth rates, one may also note 
that the range of non-linear behavior starts earlier on smaller waves 
than it does on larger ones, and this effect seems to be enhanced by 
viscosity. 

Bubble Release Frequency 

Up to this point we have dealt with a quantity, ou, which we have 
called either the "wave growth rate" to designate its functional meaning, 
or the "wave frequency" to indicate its relationship to a stable travelling 
wave . 


Now we wish to define an almost unrelated variable, the bubble 
release frequency , f b . As we can see from Figs. 42 , 43, 44 , 49 , 50, 
51, and 52, the linear portion of wave growth is only a small fraction 
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k Pa, q = 0.85x10 s watt /meter 2 , 
f b = 2 l bubble /sec, w dF =46.3 hertz, 
R 1 = 0.29 , R' =0.41. 
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Fig. 40 Growth of waves on a 1.03 mm dia wire 
heater in cyclohexanol. P = 8.05 k pa,q = 

0.98 x I0 5 watt/meter 2 , f b = 22 bubble /sec, 
a 'd/ 46 - 5her,z > R ' = 0 295, R' =0.41 
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Fig. 41 Growth of waves on a 0.406 mm dia 
wire heater in Iso - propanol. P=98.5 k 
pa,q = 1.7 x|0 5 walt/meter 2 , f b = 31 
bubble/sec, u> d = 50.0 hertz, R 1 =0.13, 

R' = 0.17. F 
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Fig. 42 Growth of waves on a 1.295 mm dia 
wire heater in cyclohexanol. P = 98.5 
kPa,q = 0.5l4x I0 5 watt/meter 2 , 
f b = 2l bu bble /sec, cu dfr = 50.0 hertz, 

R 1 = 0.425, R' c = 0.54. 
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Fig. 43 Experimental observation of frequency of 
dominant wavelength on cylinders. 




of the entire growth process. Furthermore, we haven't even attempted 
to record the recovery of the interface after the bubble leaves. It is 
clear that the periodicity of bubble growth is very weakly related to uu. 
Indeed we have already noted that Berenson's major contribution to the 
hydrodynamic theory of Q m ^ n was to reject Zuber's suggestion that f^ 
might be expressible as a simple function of a ! . 

Our wave growth rate observations give us both f b and w for 
different wavelengths and liquid viscosities. We plot this information in 
Fig. 44. The wavelength and the liquid viscosity seem to strongly 
affect the ratio of bubble release frequency and linear growth rate. 

This effect is introduced by the non-linear growth of the wave. The 
chained line shows the value of f^/w assumed by Zuber in his pre- 
diction of the minimum heat flux on flat plates. The dotted line shows 
the value of f b /w as proposed by Berenson. He fixed this ratio so that 
it was consistant with his observations of the minimum heat flux on 
flat plate heaters. 

Clearly, there remains much to be learned about the general 
dependence of f b on M and R' , and we have not even addressed the 
very important influence of varying heat flux on f b in this study. 

Effect of High Volume Fluxes on the Dominant Wavelength 

In the previous sections we saw that wavelengths longer than the 
"most susceptible" wavelength occur frequently in real systems. These 
waves obeyed the dispersion relation well, but they did not assume 
exactly the maximum frequency. Since most of our experiments were 
conducted at very low pressures, we presume that the longer, slower- 
moving waves are preferred when large vapor volume fluxes are genera- 
ted. This presumption is reasonable since we also find a stretching 
of the wavelength when the heat flux is substantially above q^^ (see, 
e.g . [63] and a discussion of these results in [5]). 

It also turns out that we might well expect to find wavelengths 
shorter than under certain conditions in which the volume fluxes are 
very small. The film boiling studies of carbon dioxide near its critical 
point by Abadzic, Grigull^and Goldstein [64, 65] show the presence of 
nearly stationary waves with wavelengths close to the critical value. 

In this case the volume flux is quite low. 

In Figure 45 we display dimensionless wavelength as a function of 
volume flux for horizontal cylinders. Wavelength data have been 
obtained both at low pressures and high heat fluxes and are given in 
Table VII. The coordinates have been reduced to dimensionless inviscid 
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Fig. 44. Effect of viscosity and wavelength on bubble 
departure frequency. 
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Fig. 45 Effect of volume flux on the dimensionless dominant 
wavelength. 





Table VII Data Used to Obtain Effect of Large Vapor 
Volume Fluxes on Dominant Wavelength 


R 

, X i 

(mm) | 

l 

R' 

i 

R' 

c 

M 

1 

\>x 10 -2 
(m^/fn^hr) 

Obse rved 
(minimum ■ 

Wavelength 
- maximum) 

X (mm) 

A 


Cyclohexanol 

j 

■ 


0.203 

0. 1 1 

0.18 

1 

16 

21 

8.89-10.67 

0.44-0.53 

1 



f 

1 

31 

10.67-13.46 

0.53-0.67 

0.514 

0.276 

0.36 


104 

16.00-21.34 

0.80-1 .06 



- 

+ 

J 

: \ 

118 

18.80-25.40 

0.94-1 .27 

: 


0.280 

0.35 

52 

3.7 

12.70-15.20 

0.67-0.80 



f ; 

i 

i 

t 

5.8 

14.73-17.27 

0.77-0.91 



i 

0.38 

16 

23 

,18.03-19.56 

0.90-0.97 

. 

■ 


} i 

l j 

0.39 

35 

7.8 

13.97-19.56 

0.71-1.00 



0.290 i 

0.41 

t 

8.8 

13.97-16.76 

0.73-0.88 

1 

* 

0.295 

0.42 

j 52 

; 5.0 

1 13.97-17.02 

i 

0.73-0.89 

! 

Acetone + 

. l . _ i * - ■ 

— 

jo. 152 

0.101 

0.13 

i 599 

J 

1 

| 0.40 

I 

2.69- 4.83 

j 

i 

0. 16-0.28 

0 . 267 

0. 161 

0.19 

j 

j 0.42 

3.88- 6.78 

0.22-0.39 

jo. 41 3 

0.257 

0.30 

j i 

! 0.62 

: 6.56-11.12 

0.38-0.64 

0.654 

0.458 

0.60 


i 0.70 

j 9.40-14.24 

0.54-0.82 

0.826 

0.514 

0.59 


| 0.97 

1 10.82-15.29 

0.63-0.89 

1.028 

0.650 

0.70 

* 

i 

j 0.82 

12.44-18.41 

0.72-1 .07 


+ Data reduced from observations of Sun as reported by Lienhard 
and Carter [1]. 
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flat plate values so that the data obtained at different situations can 
be compared. The vapor volume flux has been normalized with 
Lienhard and Wong's correction for the effect of radius of curvature 
of the cylindrical heater on the minimum heat flux. Similarly the 
wavelength data have been normalized with their expression representing 
the effect of heater radius on the wavelength. 

There is a definite increase in the dominant wavelength with 
increasing vapor volume flux. Photographs of the film boiling process 
at different heat fluxes are shown in Figure 46. The longer wavelengths 
at the larger volume flux in the bottom picture can easily be distinguished 
from the shorter wavelength in the top picture. 

On Predicting the Minimum Heat Flux 

We should now like to refer back to pages 7 and 8 in which we 
discussed Zuber's prediction of q^p for a horizontal flat plate. 

Equation (6) was the general starting point for any q min prediction. 

For the horizontal flat plate, we can accept everything that Zuber 
and Berenson did up to equation (7). Using the notation of the present 
Chapter in which tD is a real quantity — the exponential growth rate, we 
have : 


ci . = constant h x (j a \ 

min r f fg a F max 

There are two difficulties with this result. The first is that 
will increase if the fluid is very viscous. The second is that it is ^ 
based on the presumption that the minimum f fa is proportional to uu 
The previous section on bubble frequency cast some doubts on this 0 "* 3 * 
idea and it showed that viscosity tends to slow U). Thus q mi could 
either increase or decrease with viscosity. The same sort of difficul- 
ties are inherent in predictions of q^^ for cylinders. 

Furthermore, various previous authors [66], [5] have pointed out 
that many physical variables, such as material of the wire, the end 
mountings, the material and size of the wireholders, and the way the 
end mountings are attached to the wire, can seriously affect the 
observations of minimum heat flux. Data obtained under different 
conditions have been found to differ with each other by as much as 
100 to 200 percent. With so many variables affecting the true location 
of the minimum heat flux, an elaborate experimental program would be 
needed to provide reproducible and meaningful results. Therefore, we 
did not make any special effort toward measuring the minimum heat 
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Stretching of wavelengths at higher vapor volume fluxes 





flux in the present work. 


We thus leave the question of predicting q min as an open one. 
At this time it is very difficult to say quantitatively what the net 
effect of viscosity is on the minimum heat flux. 

Conclusions 


1. ) An analysis of the Taylor stability of both plane and cylindri- 
cal interfaces has been completed. 

2. ) An increase of wavelength with the liquid viscosity has been 
measured and found to compare well with the theoretical predictions. 

3. ) The wavelength and frequency measurements are faithful 
to the predicted dispersion relation. 

4. ) A successful correlation for the vapor blanket thickness 
around a cylindrical heater has been established. 

5. ) The bubble release frequency seems to depend upon both the 
dominant wavelength and the liquid viscosity. 

6. ) The generation of large volume fluxes can cause the wave- 
lengths in film boiling to considerably exceed the "most susceptible" 
wavelength. An experimental correlation of wavelength with volume 
flux is presented. 

7. ) The minimum heat flux has been shown to be sensitive to 
licfjid viscosity. The nature of the influence of viscosity on q . is 
not presently known, however. 



V. THE PEAK HEAT FLUX IN VISCOUS LIQUIDS 


The peak heat flux on any heater configuration is given in equation 
(3) as 


q = p h (A ./A.) U 
max K g fg v j h' 


H 


we shall for the moment restrict our attention to the flat plate configu 
ration and ask how Aj/A h and U H change in viscous liquids. 

We saw in the last Chapter that the wavelength will stretch with 
increasing viscosity. However, the ratio (Aj/A^)^ is independent of X^ 

and will equal tt/ 16 for any flat plate. The ratio A^/A h may vary with 
X d in certain other configurations but that is a matter we shall defer 
for the moment. 

The important variable in the problem, and the more difficult to 
predict, is U^. We shall attempt first to estimate U H for a jet leaving 
a flat plate. 

The Maximum Stable Vapor Velocity in a Jet 

Velocity Profile for the Primary Flow in the Vapor Jet and in the 
Liquid Column . Figure 48 contrasts the vapor removal idealization for 
the inviscid case with that which we shall assume for the viscous case. 
In the viscous case we must avoid the configuration of a circular jet 
in a square area if the analysis is to be at all tractable. Thus we 
inscribe a circular liquid return path inside the square area. We also 
make the following assumptions: 

1. ) The gas velocity, U g , and the liquid velocity, U f , are zero 
on the interface. 

2. ) The PoiSeuille velocity profile is valid in the vapor jet. 

3. ) The analysis is two-dimensional. However curvature of the 
jet will be introduced subsequently as an additional pressure contribu- 
tion and velocities will be adjusted such that: 

4. ) the mass inflow of liquid in the circular annulus is matched 

by the gas outflow in the circular jet. U^ is locally maximum or 
minimum at y = 2Rj. 

5. ) The shear stresses are matched at the interface. 
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The resulting velocity profiles are obtained by fitting even poly- 
nomial functions in y/R. to these conditions. For the gas we obtain 


U 
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U 


(1 - 


m 


y 2 /^ 2 ) 


(89) 


where the subscript m denotes the maximum value; and for the liquid 
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(90) 


These results incorporate considerable physical approximation and, 
indeed, when we complete the analysis we shall find it necessary to 
adjust one constant before the theory will fit the experimental results. 
However, the actual physical configuration is so complicated that we 
could hope for no more. What the present analysis will yield is the 
proper form of the dependence of the peak heat flux on the relevent 
dimensionless groups. In the light of this objective our failure to 
represent the configuration more precisely is of no great importance. 

Equation (90) gives a maximum downward velocity of: 

p o 

+ 0.0388 j (91) 

Pf " 


U. 


= -u. 


m 


m 


0. 


3557 


M 

IS 


which occurs at y = 1.28 R.. We shall need these results subsequently 
in the analysis. It is wortp} noting that this maximum occurs near the 
interface and that there is a slight upward recirculation of liquid at 
y = 2R.. We should also be aware that the constants in equation (91) 
fully reflect any crudity in our assumptions as to the liquid and vapor 
flow configurations. 

Formulation of the Stability Problem. To study the stability of 

the interface we superpose perturbation components of velocity and 

pressure on the mean flow of gas. We nondimensionalize velocity, 

distance, pressure, and time with U F , R., p_ U f 2 , and Ri/U f , 

r m J 9 m J m 
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respectively. We then define, u, v, and p as the dimensionless per- 
turbation components of axial velocity, radial velocity, and pressure. 

And we use the subscripts t, x, and y to denote differentiation with 
respect to the dimensionless dependent variables. 

The continuity and Navie r-Stokes equations for the perturbed flow 
on the gas side are then: 

u + v =0 
x y 

u + Uu + U v = -p + Re"\u + u ) 

t x y K x v xx yy' 

v f Uv = -p + Re Vv + v ) 
t x y xx yy 

where Re = U frn R P At . The boundary conditions at the interface are: 

J y 


(92) 

(93) 

(94) 


v = 0 

at y = 0 (95) 

u =0 

y 

p - 2Re v + cm. /(p LL 2 R.) + art '(2p LL 2 R.) 

9 y xx 7 g fm y x v g fm y 

at y = 0 (96) 


= P F at y = 1 (97) 

' complex 

u v/ + v = Re r f at y = 1 (98) 

y * complex 

where ti Is the displacement of the interface between the gas jet and the 
liquid column from its mean position. p f and are 

complex complex 

the complex perturbation pressure and shear stress respectively from 
the liquid side. 

Equations (92) through (94), with boundary conditions (95) through 
(98), have been solved for small values of Re and for wave numbers 
of the disturbance close to 1. Complete details of the analysis are 
given ip reference [27]. For a neutral disturbance with a wavelength 
equal to the "most susceptible" Taylor wavelength in the interface of 
the gas jet and the liquid column, it is shown in reference [27] that:^ 


q 

The solution is fairly complicated. It involves putting the Naviei — 
Stokes equation into the Orr-Sommerfeld form and obtaining the lead 
terms of a series solution to it. 
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1.07 


CT 


(99) 


II -ft 

11 11 2Rj 

where represents the magnitude of the imaginary component of the 
complex fluctuating shear stress and p^ represents the real component 
of the complex fluctuating normal stress that results from the down- 
ward flow of the liquid. 

The evaluations of t ^/r\ and p^A) are extremely difficult. They 
involve solving the same governing equations with appropriate boundary 
conditions. Brooke Benjamin [67] has solved the comparable problem 
for a boundary layer type of exterior flow on a way wall, with no 
pressure gradient. For the shear stress and pressure perturbations, 
he obtained: 


2 r 2 

— R; P f U f 

U 1 f f m 


( 100 ) 


r f x 2HR. v.5/3 

— = -1.188' ^ ) p 


f Uf m 2 ( 


Mf 




1/3 


T u f> R j 2 >y=R j 


( 101 ) 


We shall use these expressions in equation (99) with the realization 
that they represent the incorrect configuration. The liquid velocity 
profile in our case is more diffuse and it has a negative pressure gra- 
dient in the direction of flow. However, if we did obtain equations 
equivalent to (100) and (101) for our situation, they would differ only 
in the lead constant on the right-hand side. Substituting equations 
(100) and (101) in equation (99) and rearranging the result, we obtain: 


U 


2 

f l 1 

m 


- 1/3 


-1 .63 


L p f (Uf) R 2 
f f y j 




TTRjPf 


( 102 ) 


y=R j 


To get the desired expression for the maximum stable vapor veloci- 
ty in the jet we combine equations (102) and (91) and replace R. with 
A^/4. The result is ■' 
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f d 


The form of equation (103) is correct within the limitations of the 
three constants 0.3557, 0.0388, and 1.01, which reflect the crudities 
of the various configurational assumptions. 

Peak Heat Flux Prediction 


The average vapor velocity for Poiseuille flow in a cylindrical jet 

at the point of instability is U /2. Using U n /2 for U u in equation 

m om ’’ 

(3) and A^/ = tt/ 16 we obtain: 


max. 
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Substituting equation (103) for U with p /p ~ 0, and dividing by q 


1.14 q 




maxr 


max. 


we obtain the desired expression for the dimensionless 

viscous peak heat flux on a flat plate. (This assumes p /p « u /u ' »• 

g f g 'V • 


max p 


(q ) 

r mviscid 
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Our approximations as to the flow configuration are now reflected 
entirely in the two constants: 0.633 and 0.807. 


In addition to the now— familiar dimensionless groups, 
we have now added a new one: 


M and A^, 


V = 


v' Pg Pf a 


3/4 


t* tg(Pf-p g ) 3 
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It is very interesting to note at this point that Borishanskii [36] 
suggested in 1956 that a variety of up to 8 dimensionless groups might 
be involved in the viscous peak heat flux problem . These included 
groups which could be combined to give both M and V. However, he 
suggested that a group equivalent to M, alone, would be sufficient to 
characterize the viscous effects. 


In 1970, Lienhard and Keeling [4] identified two types of viscous 

influences on q One of these was the result of induced convection 

m sx 

which we have mentioned briefly on pages 12 and 13, and which we 

carefully avoided in the present experiments. The second is the 

direct influence of viscosity on the components of the hydrodynamic 

processes. Both Lienhard and Keeling, and Borishanskii sought to 

correlate data in terms of M (and R' in the case of finite 

max Ml dX —7 


geometries) without using additional groups. They both succeeded, but 
in each case they were restricted to values of M significantly greater 
than 100. Borishanskii's equation (10) was simply based on insufficient 
data, and values of his N (approximately equal to our M 2 ) that were 
too large to reveal viscous influences with any certainty. 


We have evaluated equation (105) numerically and plotted the re- 
sults in Fig. 49 for four typical values of V, with M as the independent 
variable. In order to do this, we have used the R' = “curve from 
Fig. 24 to express as a function of M. By following the saturated 
liquid-vapor line for one of our experimental fluids — cyclohexanol — at 
earth normal gravity, we have obtained a unique value of V for each 
value of M. Thus Fig. 49 displays a locus of points for cyclohexanol 
on dimensionless heat flux vs_ M coordinates. This particular curve is 
similar in form to equation (10) although it is higher. 


Comparison of Theory and Experiment 

We have already described experiments in Chapter II which genera- 
ted new data for flat plates, free of induced convection effects (recall 
Table I on page 20). Table VIII presents additional data obtained in 
this apparatus using cyclohexanol under low pressure. These data 
embrace values of M, sufficiently low to reflect a significant influence 
of viscosity. 


Since we have already seen in Fig. 49 that, for a particular fluid, 


at one value of g, equation (105) takes the form q 


'max, 


/(q, 


max. 


) 


V r— III • 

t" r inviscid 

= f(M, only), we plot the viscous cyclohexanol data for earth-normal 
gravity on these coordinates in Fig. 50. These data have been corrected 
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Table VIII. Peak Heat Flux Data for Cyclohexanol 


for Flat Plate Heaters 


9/9 e 

P 

— 

M 

! v 

“ | ? 

L* 

q x10 5 

q 

max r 

F 

(kPa) 

t 2 tt/3 

! 

max 

(W/m 2 ) 

(q ) 

v maxp'' 

inviscid 

1 

1 .61 

20 

260 

3.23 

1 .67 

2. 19 


3.62 

35 

337 

3.31 

1 .99 

1.75 


5,66 

48 

380 

3.35 

2. 10 

1.54 


11.73 

100 

490 

3.45 

2.55 

1.34 

V 

22.52 

155 

600 

3.55 

2.81 

1.25 

7 

i 

99.50 

, 250 

920 

I 

| 

3.86 

5.65 

1 .02 

4.97 

11.86 

I 70 

1 437 

7.69 

3.75 

1 .48 

8.95 

7.31 

38 

i 

350 

J 

10. 1 1 

_J 

3.53 
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Fig. 50 Effect of liquid viscosity on peak h^at flux for both 
flat plates and cylinders in cycloh.xanol at earth 
normal gravity. 



for the effect of finite plate size since fewer than 3 wavelengths are 
accomodated on the plate (recall Fig. 8) after the wavelength is cor- 
rected for viscosity (recall Fig. 24). 

With these data we include original earth-normal gravity data for 
cyclohexanol boiling on horizontal cylindrical wires of various sizes and 
at various pressures. The data have been nondimensionalized by divi- 
ding them by the corresponding inviscid values. These data were ob- 
tained using exactly the same techniques as those employed by Sun and 
Lienhard [6], and they are presented in Table IX. 

The comparison of data for cylinders with data for flat plates is 
justifiable if R' is near unity, and we use data in the range of 0.22 s 
R' ^ 1.24 here. Sun and Lienhard's [6] inviscid hydrodynamic model 
for the peak heat flux on small cylinders assumes the gas jets are 
separated by A. d . They took the jet radius to be (R + 6), where 6 is the 
thickness of the vapor blanket surrounding the heater near the peak 
heat flux. While the radius of the jet decreases as (R + 6) ? so too does 
for small cylinders. The effect of a decrease in the liquid return 
velocity due to the smaller jet radius is nearly compensated by this 
shorter Thus this model differs from our assumed flat plate model 

only in that the location of jets is on a line rather than on a square grid. 

Figure 50 also includes equation (105) for comparison. However, 
as we noted in deriving this expression, the two decimal fractions in 
it should be regarded as free constants which might have to be altered 
empirically. Of these two, the value 0.807 in the denominator is 
relatively inconsequential, since it multiplies a term that is generally 
fairly small. Thus we shall not tamper with it. On the basis of the 
data in Fig. 50, however, we shall change the value 0.633 in the 
numerator to 0.264. With this adjustment, the comparison between 
theory and experiment in the figure is very good indeed. 

When we go to variable gravity experiments in cyclohexanol, V 
can vary independently of M. Accordingly ‘Vnax^Snax^inviscid becomes 
a function of two variables, M and V. Equation (105), with the lead 
constant equal to 0.264, has been plotted as a three-dimensional sur- 
face in Fig. 51. Fifteen high-gravity data for horizontal cylinders 
(Table IX) and one high-gravity data point from the flat plate apparatus 
(Table VIII) are included in the figure. 

We see that these data faithfully bear out the independent influence 
of the parameter, V, as it is predicted by equation (105). 


Table IX. Peak Heat Flux Data for Cyclohexanol Boiling 
on Horizontal Cylindrical Heaters 


R 

(mm) 

R' 

9/9 
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P 

(k Pc.) 

M 

V 

ci x 1 0 0 

vriax 

( 5 ) 

^ r> ret 

(9 ) 

max , . . 

inviscid 

0.41 
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47 
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61 
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4 • 26 
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14 
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3.25 

4.34 
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2.76 

30 
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3.37 
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67 
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3.69 

2.01 
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8 

,5 

3.54 

20 
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18 
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It is important to note that as M increases or as V decreases, 

the prediction of q~ av /(q~*=~) eventually falls below unity. 

max max inviscid 

At this point we simply abandon the viscous prediction and embrace 
the inviscid prediction. This is reasonable since a basic limitation of 
any viscous theory is that it does not pass to the limit of inviscid 
behavior as (I is made very small. 

Conclusions 


1.) We have predicted that 


q 

maxp 

Constant 

V/MA d 1/2 

(q ) 

v maxr: . 

r inviscid 

y 1 -0. 807 r V 

inviscid -V3 


y Lm 2 a 

q J 


“VI ‘ , 

d 

niaxp 

for viscous liquids. 




2.) The Constant has been evaluated experimentally. It is ~ 
0.264. 


3.) This prediction can also be applied to data for cylinders as 
long as R' is on the order of magnitude of unity. In this case, 

/(q ) should be rewritten with the appropriate values 


max. 


max 


F inviscid 


of and (q ) 

max max inviscid 


for cylinders. 


4.) The viscous prediction is only valid when it predicts higher 
peak heat fluxes than the equivalent inviscid theory. 


5.) The derivation of equation (105) required the solution of an 
original stability problem, namely the prediction of the maximum stable 
velocity of a gas jet in a viscous liquid. This velocity is given by 
equation (103). This problem might well find applications beyond the 
present one . 
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VI. HYDRODYNAMIC THEORY OF ELECTROLYSIS 


Introduction 


This chapter will deal with a non-boiling application of the prin- 
ciples of the hydrodynamic theory. We shall discuss the gas removal 
process during high current-density electrolysis with the aim of 
identifying hydrodynamic transitions, analogous to those which occur 
in boiling. This work was extracted from the Master's thesis of 
Bhattacharya as reported in [18] and [19]. 

Electrolysis has attracted increasing interest lately because of 
its relevance to a variety of problems of human life support and a 
myriad of industrial processes. References [68] and [69] provide some 
background as to the scope of contemporary applications of electrolysis. 
This particular study was part of a broader effort to establish the 
similarity between boiling and electrolysis. Bhattacharya generally 
sought to establish means for applying known methods for predicting 
boiling processes to the prediction of comparable processes in 
electrolysis. An earlier paper from his study, reference [17], dealt 
with bubble growth predictions for the two processes. 

Figure t>2 is a conventional boiling curve with a transformed ordi- 
nate. In it, the vapor volume flux, v, (or the heat flux, q, divided 
by the volumetric latent heat of vaporization, p h f ) is plotted against 
the driving temperature difference. At. The \a^ 11 fmown transitions 
in the boiling mechanisms are indicated on this figure. Let us consider 
each one briefly: 

1. ) inception: This is the point at which the first bubbles appear. 

Since we are not concerned with predicting inception in either boiling 

or electrolysis, no attempt will be made here to recount the massive 
literature on the subject. 

2 . ) transition from the isolated bubble region to the region of 

slugs and columns: This transition was identified by Zuber([70] and 

[71]). It was subsequently predicted by Moissis and Berenson [32], who 
also advanced a variety of measurements of the transition point. 

Using observations of the rise of bubbles in tubes [72] and [73]. 

Moissis and Berenson shewed that the volume flux at this transition, 

v , can be given as 
sc M 
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observed during boiling 
conjectured for electrolysis 
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0 driving temperature difference, AT, or voltage, E 

Fig. 52 Schematic boiling curve and corresponding conjectured electrolysis 
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V 

sc 


0.56 


3 1/2 4 /g ~/(P ~-P ) 

f g 



(107) 


where 3 is the contact angle of the liquid— vapor interface on the 
heater. Zuber's value of A./A^ = tt/ 16 for flat plates in this expres- 
sion led to 


v 

sc 


= o.n 3 


1/2 4 


P f“ ^g) 


(108) 


In the published discussion of [32], both Chang and Zuber asked 
the authors if this transition was well defined. They replied that it 
was indeed slurred so that equation (108) is an approximate result. 

We shall return to this point later. 

3 *) the peak volume flux transition : Dividing Zuber's and Sun's 

expressions for q max by P g hfg we obtain, respectively: 

■ tt 4 ! — — - 2* 

V max ~ ^ v'9° ^’f~ P g)/ p g (109) 

and 

V max , ~ v max (9*39+2.27 exp [-3.44R* ]) (110) 


4 0 the minimum heat flux transition : We likewise obtain from 

Zuber's and Lienhard-Wong's [33] q m i n expressions, respectively: 


and 


v . = C 

min 1 



°9 


(0.+ o ) 
f g 


V 

min 


cyl 


v 

mm 




_R 7 2 (2R /2 + 1) J 


1/4 


( 111 ) 


(1 12 ) 


where C-| was predicted as 0.171 and C 0 was found experimentally to 
be 1 . 289 . 
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There is a general failing to these Vp^p expressions that does 
not afflict As we noted earlier, Berenson [35] showed that 

while v max is insensitive to surface condition, v min can be strongly 
influenced by minor variations of surface condition. Kovalev [36] 
further amplified the point and showed that the end-mounting of cylinders 

also influenced v m ^ n . On the basis of Kolavev's data, Lienhard [7] 

cyl 

showed that C 2 should be reduced to at least 0.0217 — perhaps even 
further. Thus v predictions are generally a chancey business. 

What the real minimum is, cannot finally be predicted, and it is very 
hard to ascertain experimentally. 

The Transitions in Electrolysis 

At first glance we should expect the same transitions to exist in 
the volume flux vs. voltage curve for electrolysis, as exist in the 
boiling curve (see Fig. 52). While v is equal to q/Pghf boiling, it 
can be obtained from Faraday's law for electrolysis. Faraday's law 
expresses the mass flux of Hg dissociated at the cathode as proportional 
to the current flux, 1/A amps per unit area, the molecular weight (M) of 
H 2 , and n -1 where n is the number of electrons transferred in the 
liberation of one H 2 molecule (for H 2 , n and M both equal 2.) Thus: 


v = (I/A) M 
96500 o n 
H 2 


cm 3 


2 

cm 


sec 


(113) 


and we can reduce the current to a volume flux just as we could q. 


Let us now consider a basic difference between nucleate boiling 
and nucleate electrolysis. The basic mechanism of heat removal in 
nucleate boiling was shown by Zuber [70] and [71] to be microconvection 
driven by the rising bubbles. In 1962 Tien [74] predicted the micro- 
convective heat flux in the region of isolated bubbles using an inverted 
stagnation flow model. Boehm and Lienhard [ 75 ] summarized and 
extended this type of prediction in 1964 and found that in the region 
of isolated bubbles 


.1.18 1/3 

q ~ v ~ A T n 3 (114) 

where the "site density", n s , is the number of nucleation sites per unit 
area. 


During nucleate boiling the site density typically rises as 
Ay 5 or more . Thus q or ^ typically rises as At 3 or 4 . But during 
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electrolysis there should be no corresponding convective enhancement 
of charge removal. The diffusion coefficient for H created at the 
wall is very small (about 1/30 of the thermal diffusivity of the water). 
Thus the Schmidt number is on the order of 100 and the diffusion 
boundary layer is very much thinner than the microconvective velocity 
boundary layer. Furthermore, we find a heavy density of idle bubbles 
scattered in among the active or repeating sites during electrolysis. 
These create a static sublayer of fluid whose thickness is independent 
of the microconvective velocity, and which should totally contain the 
diffusion layer. Thus v should simply increase in direct proportion 
to the driving voltage, E. None of the convective reinforcement by 
the nucleation sites, which occurs in boiling, is present. In electroly- 
sis it is as though the sites were not there, insofar as the relation 
between v and E is concerned. 

An important aspect of the transition from the region of isolated 
bubbles to the region of slugs and columns is the blanketing-over, 
and removal of liquid-metal contact that occurs. Throughout the 
isolated bubble region there is virtually no loss of contact; but as 
the slugs and columns evolve, there is. Zuber [71], in fact, envisioned 
that this loss of contact increases, honeycombing the liquid under the 
jets and feeding them. This can continue until more than tt/16 of the 
contact is eliminated and large dry patches cover the heater (see, e.g., 
the photographs of Kirby and Westwater [76] which dramatically show 
this behavior). 

The heat transfer continues to increase strongly with At in the 
region of slugs and columns because as the vapor velocity increases and 
heater surface is exposed, the microconvection at the surface continues 
to improve. Quite the opposite should occur during electrolysis, 
however. As contact is lost, the electric current — and with it, v — 
must decrease in the absence of a convective process. 

We are interested, then, in learning to predict the hydrodynamic 
transition from the region of isolated bubbles to the region of slugs and 
columns, or the "first transition" as we shall call it. This transition 
should replace the Helmholtz unstable transition as the point of peak 
volume flux, in electrolysis, it can be predicted for electrolysis on a 
flat plate with the help of equation (108). Since the experimental work 
presented here is done on horizontal wires, it will next be necessary to 
develop an equivalent expression for cylindrical cathodes: 

On a cylindrical element, Aj/Ap, must be modified for use with 
equation (107). The bubbles rise from the wire, as indicated in Fig. 53 , 
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vapor jets-. 
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eventually following one another without interruption in a column of 
diameter, D^, equal to that of the departing bubbles^. In this case 

Aj_ = (n /4) D b 2 _ D b 

A h 2TTRL 8C 3 R C 


where L, the length of wire subtended by one column can no longer 
be assumed equal to X- ., since we are far below \> = . It is therefore 

replaced with CgD^, where Cg is an unknown constant. The constant 
1/Cg can be viewed as the fraction of the wire's length occupied by 
columns. In Moissis and Berenson's model, Cg would be equal to 2, 
measured along the rectangular grid lines, and equal to 2/2 or 2.83 
on diagonal lines through the grid. The spacing should be comparable 
for cylinders, but we have no way of specifying the precise value of 
C 3 a_ priori . 

Using Fritz's expression (quoted in [79]) for the bubble departure 
diameter, D b , 

! o n 

D h = 0.0148 3 — (116) 

b 9< P f" P g) 


in equation (115), and then substituting equation ( 115 ) in equation ( 107 ) 
we obtain 


v 

sc 


0.001464 



g3/2 

R' 



(1 17) 


for the volume flux at the first transition on cylinders. 

Ex periments 

A local maximum in electrolytic gas evolution has been observed 
in the past. In 1950, for example, Kellog [ 78 ] reported observations 
of such a maximum during aqueous electrolysis from a vertical wire 
anode. Actually this transition was known then in connection with the 
"anode-effect", or the vapor blanketing of anodes which had been 
observed much earlier during the electrolysis of molten salts. The 


10 

This conceptualization is consistent with that of Moissis and 
Berenson for flat plates. 
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concern in this case was not with predicting (or even explaining) the 
peak but with identifying the phenomenon of film electrolysis and, 
noting that it resulted in a considerable reduction of current at a 
given voltage. One important en passent observation made by Kellogg 
was that "the current increases almost linearly with increased voltage" 
in the region approaching the local maximum . This concurs with the 
result we anticipated in the previous section. 

Two basic experiments have been undertaken in the present study 
to test our predictions as to the nature of the electrolysis curve, in 
one of the tests we measured the 1 (or v) vs . E curve during electro- 
lysis from cylinders, and obtained photographic verification of aspects 
of the process. In the other we measured v sc during boiling to provide 
additional corroboration of equation (117). 

Figure 54 shows the apparatus for the electrolysis experiment. 

A one-tenth normal KOH electrolyte was contained in a 9 cm x 10 cm 
x 19 cm plexiglas box. Clean nichrome wire cathodes (14, 16, and 20 
gauge and about 5 cm in length) were suspended in the electrolyte. The 
anode was made of clean platinum foil and suspended in the electrolyte 
above the cylinder. It was made considerably larger in surface area 
than the cathode to assure that it would function in the isolated bubble 
regime long after the cathode underwent any transition. Additional 
details relating to the experiment are given by Bhattacharya [11]. 

Figure 55 shows the basic output of the experiment for three of 
the wire sizes (14, 16, and 20 gauge). These v vs . E curves were 
obtained by increasing the voltage across the test cell up to the desired 
value, stabilizing the temperature within an 8°C range, and measuring 
the current, l. Equation (113) was then used to calculate v from l. 

It was necessary to operate in a fairly high temperature range because, 
at high currents, the capsule heated up rapidly during the course of an 
observation. We were generally able to hold temperatures close to 
70°C when the potential was below about 100 volts; but it was increas- 
ingly hard to do so when we had to increase the voltage through the 
region of peak current before making an observation. 

The curves were quite regular and reproducible up to a maximum 
value of v. Beyond this maximum, the data exhibited a great deal of 
scatter and indeed the current was unsteady and tended to fluctuate at 
a given voltage. Nevertheless the current clearly dropped off drama- 
tically with voltage beyond the peak. At the maximum volume flux a 
blue corona first formed in the vapor separating the cathode from the 
liquid. As the voltage was further increased the corona also increased. 
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Fig. 54 A test capsule and the circuit diagram for the electrolysis experiment 
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during electrolysis. 



and with it there developed intermittent sparking between the cathode 
and the electrolyte. As we approached our maximum voltage capability 
of about 180 volts, the sparking assumed a regular spacing and 
fairly great intensity. 

We should emphasize that the abscissa of Fig. 55 (the voltage 
between the cathode and the anode) is only proportional to — not equal 
to — the driving voltage of the surface electrolysis process. This 
fact was clarified by a duplicate run on a 16 gage wire using a 
reference electrode located 13 mm below the cathode. The reference 
electrode reveals that about 75% of the voltage drop is being sustained 
across the cell, away from the cathode. The actual surface voltage 
drop would be still substantially less than the value, E c _ b , measured 
on the reference electrode. 

Our system was therefore one in which' the current varied almost 
entirely by virtue of changes in potential across ir.c cathode, the 
remaining resistance being constant. The curves accordingly should 
have the right shape although the abscissas are arbitrary within a 
constant factor. 

An additional source of error in the system is the possible contri- 
bution of vapor generation due to heating. Since this could hardly have 
been a significant factor before the onset of sparking and corona, and 
since sparking only occured once the peak had been reached, its contri- 
bution to our observations of the peak could not have been great. 

However a second source of error is compensatory and comparable in 
magnitude. This is the fact that "current yield" efficiencies can run 
below 100% (see e . g . [79]) so that Faraday's Law might underestimate 
the hydrogen gas yield by a few percent. These effects probably give 
a net error somewhere between 0 and 10%. At voltages beyond the peak 
there might be considerable vapor generation. Probably the film 
electrolysis regime is, in fact, a mixed boiling and electrolysis region 
owing to the greatly increased surface resistance at the cathode and 
the resulting heat generation. 

The results of the tests as plotted in Figs. 55 and 56 verify the 
essential linearity of current (or v) in voltage, after nucleate electrolysis 
becomes fully established. This linearity is evident only after a layer 
of bubbles has been established to preclude microconvection. 

Figure 57 is a sequence of pictures of electrolysis on the 16 gauge 
wire (R' = 0.278) beginning slightly below the maximum current and 
continuing to E = 176 volts. The 1.1 cm wide marker is evident at 
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v (cm 3 /cm 2 hr) 



Fig. 56 Comparison of 1(E) and I(E c _ b ) on 
logarithmic coordinates. 





108 volts 






the bottom of each picture. The corona and sparking are invisible in 
these pictures, back-illuminated with a high -intensity stroboscope. The 
slugs and columns are very clear on the wire in the first two pictures. 
As E is increased these jets appear to give way to a film process 
characterized by Taylor unstable waves very similar to those observed 
during pure film boiling (see, e.g . , [ 33 ] or [ 5 ]). 

The second experiment used Sun's [ 6 ] apparatus to observe q at 
the first transition on horizontal wires during boiling in methanol and 
water. Since the apparatus and basic procedure do not differ from that 
described in [ 6 ] we shall not reproduce details here. In this case 
the first transition was observed with the help of a stroboscope and 
checked our visual impressions with still photographs. 

Table X lists the observed transition points for the two liquids 
and for the four wire sizes. Bhattacharya was generally able to 
identify the transitions reproducibly within a 6 % margin, but only after 
identifying a particular stage (within a broader transition range) in 
which the first few jets appeared to take form. 

Discussion of Results 


We now have the following claims as to the character of electro- 
lysis, and data against which to test them: 

1 •) Microconvection does not abet the current (or volume) flux 
in electrolysis . A corollary of this claim was that, in the nucleate 
electrolysis regime, I or v should be approximately linear in E (or 
E c-b) U P to the point of hydrodynamic transition. Figures 55 and 56 
clearly bear this out very well, as did Kellogg's experiment. Given 
that this is true then we expect that: 

2 0 The volume flux maximizes at the first transition, instead 
of the second, in electrolysis. Four kinds of evidence bear this out. 
The first is the lack of microconvective enhancement of I. The second 
is the visual appearance of the process (see Fig. 57 ) which matches the 
appearance of the first transition in boiling. 

The third piece of evidence arises from the fact that equation 
(110) predicts v max in the range from 0.975 x 10 6 to 1.19 x 10 6 
cm 3 /cm 2 hr for our test wires. Thus v max is predicted to be between 
270 to 520 times the observed maximum, depending upon the wire size. 
It follows that the maximum could not possibly reflect the second 
transition. The reader who is acquainted with boiling might justly 
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Table X 


Observed Heat Flux at First Transition in Boiling 
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wonder how this ratio could reasonably be so high. The reason is 
that for either cylinders or flat plates the ratio of transition volume 
fluxes is 


4 

V 

sc 


V 

max 



( 118 ) 


Since p g is extremely small for H 2 the ratio is quite small. The 
same would be true for boiling at low pressure. 

The fourth piece of evidence lies in our third claim, namely: 

3*) Equation (117) describes the first transition in both boiling 
and electrolysis. To apply equation (117) one must first ascertain 
the contact angle, P, for water, methanol, and tenth-normal KOH 
electrolyte, since v Sc ~ P 3 / 2 . Unfortunately P is a notoriously 
stochastic property. In 1959 Griffith and Wallis [80] pointed out the 
difficulty of reproducing P consistently. Table XI shows the great 
variability of their data with slight changes of surface condition. 
Trivedi and Funk [81] provided in situ measurements of contact angles 
under static bubbles during electrolysis in normal KOH electrolyte on 
platinum plates. Trivedi [82] recently measured P during electrolysis 
in tenth-normal KOH on nichrome wires. The data of [80] and [81] 
are included in Table XI. Both sets of data were obtained on metals 
that had been first cleaned in distilled water, then in methanol, and 
finally rinsed in the test liquid. Chang and Snyder [83] also reported 
values for methanol and water which they inferred from the data of 
earlier papers. These values are included in Table XI. 

Table XI also includes original 8 values that we measured using 
the tilting plate method ^ for methanol and the sessile drop method 
for water. Both tests used nichrome plates cleaned in the same way 
that Trivedi's surfaces were. 

The four measured values of v sc in electrolysis (the peaks in 
Figs. 55 and 56) and the eight values of v sc for boiling are given in 
Fig. 58 . They are nondimensionalized as v sc /P 3//2 Vg~d( Pf-pg) and 
plotted against 1/R' , since equation (117) indicates that this group 
should equal 0.0001 464/C3R’ • The observed values of v sc are spread 


A detailed explanation of such tests is given in Reference [80]. 
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Table XI. Measured Contact Angles 


Liauid and Metal 

Source ! 

Temp. 

C°G] 

3 range 
( degrees') 

a - = |— 

1 

water-stainless steel 

[80] 

30 

38-78 

watei — stainless steel 

[80] 

100 

28-82 

— N, KOH - platinum* 
10 

[81] 

23 

15-35 

- — N, KOH - nichrome wire* 
10 

[82] 

23 

10-23 

water - copper 

[83] 

180 

60 

methanol - copper 

[83] 

64 

46 

water - nichrome 

present 

27 

60-81 

methanol - nichrome 

present 

27 

40-52 


* data for electrolytic Hg bubbles _in_ situ . 
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Fig. 58 Comparison of predicted first- transition 
data for both boiling and electrolysis. 





out considerably in this representation owing to our ignorance of 
precise values of 3. For the boiling results we used our own measure- 
ments of 3 and for the electrolysis data we used Trivedi's values of 3 
for electrolysis on nichrome wires. 

Figure 58 shows that the 1/R' dependence of v gc is upheld, 
although the slopes differ between water and methanol. The water 
data for both electrolysis and boiling are fairly well represented with 
equation (117) if C 3 is taken to be 4.4. This value of C 3 corresponds 
nicely with the proportioning of the jets that we observe in the first 
two pictures in Fig. 57 (taken near the point of macimum v). A value 
of C 3 more nearly equal to 2.5 would be required to represent the 
methanol data in Fig. 58. 

4 .) As E is increased beyond the point of peak volume flux, 
the process goes to a kind of film process similar to film boiling . 

The photographs in Fig. 57 show an evolution from gas removal by slugs- 
and-columns to gas-vapor removal by the cyclic collapse of Taylor- 
unstable waves. The wavelengths measured from the last three photo- 
graphs in Fig. 57 and from other pictures not included here are: 

(X ) = 1.07 ± .25 cm 

° experimental 

while Lienhard and Wong’s theory predicts X d = 1.04 cm. 

The wire appears to become almost completely blanketed in the 
pictures for high voltages. However a small number of departing 
nucleation bubbles in the pictures suggests that there is still some 
liquid-to-cathode contact. Despite these indications of a conventional 
approach to a film gas removal process, there is a difficulty in the 
numerical value of v mip . Equation (112), which (as we have noted) 
might overestimate v^^, gives: 

3 2 

V fR' = 0.174) = 4550 cm /cm hr 

min v 

v (R* = 0.348) = 3020 cm 3 /cm 2 hr 

min 

while the corresponding experimental values of v sc are 3660 and 2280, 
respectively. 

Thus v sc < v m i n , and beyond E(v sc ) the calculated hydrogen gas 
volume flux drops still farther below This is additional evidence 
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that stable film electrolysis involves a substantial component of 
evaporated water which does not appear in the v calculation based on 
Faraday's Law. 30 

Conclusions 


1. ) There is an enhancement of the heat flux which occurs in 
boiling when nucleation sites influence one another; this does not occur 
in the volume flux in electrolysis. Consequently v ~ E in the region 
below the first transition. 

2. ) The volume flux reaches a local maximum at the "first 
transition" point, or the transition from isolated bubbles to slugs-and- 
columns. 

3. ) The volume flux, v Sc , at the first transition on horizontal 
cylinders is given approximately by equation (117). However, ignorance 
of precise values of 3 and of the precise structure of the escaping gas 
jets precludes accuracy beyond a factor of two in the use of equation 
(117). The predicted dependence of \> sc on 1/R' appears to be sound. 

4«) A form of electrolysis that shows many outward appearances 
of film boiling, including the Taylor unstable wavelength, are established 
at high voltage. 

5. ) A clear understanding of transitional and film electrolysis 
and an identification of v^^ in electrolysis must await experimentation 
in a system with an extremely high voltage capability and a substantial 
cell cooling system. 

6. ) The existing body of understanding of boiling processes pro- 
vides a very useful point of departure for seeking to extend our under- 
standing of gas— forming electrolysis. 


VII. BOILING ON SMALL WIRES 


Introduction 


We now turn our attention to a matter that actually lies outside 
of the hydrodynamic theory, as it has been described here. This 
chapter will describe the phenomena that replace the hydrodynamic 
mechanisms when a finite heater becomes very small or gravity becomes 
very low. These processes are extremely important in the context of 
aerospace applications since they will inevitably be involved as any large 
system passes from earth-normal gravity to the zero-gravity limit. 

Under conditions of very low gravity, or for very small heaters, 

R' (which characterizes the ratio of buoyant forces to surface tension 
forces) becomes small. Figure 59 shows the same Qppax correlation 
that we showed in Fig. 13 for Sun's horizontal cylinder data, but now 
we have plotted the abscissa on logarithmic coordinates to illuminate 
small R' behavior. It is clear that R' must exceed 0.15 if the wave 
mechanisms which define q max are to stay intact. For smaller values 
of R' , surface tension so over-balances inertia that these mechanisms 
deteriorate, and a sampling of data by Kutateladze et al. [84], Sun [6], 
and Siegel and Howell [85] no longer correlate on Pmax^max — * 
coordinates. ^ 


It was likewise shown in [5] that there is still good visual evidence 
of the wave stability mechanisms for all R' - 0.12, during film boiling 
on horizontal cylinders. But for R' < 0.06 these mechanisms cease to 
be identifiable. 


If a correlation of the form 


q m ax^max- 


= f(R') fails. 


and the 


Zubei — Kutateladze mechanisms fail with it, what exactly happens to the 
boiling process at small R'? The photographs of Sun [6], of Pitts and 
Leppert [86], and of Kutateladze et al. [84] for small wires provide 
some clue. Bubbles grow on the wire until they are large enough to 
buoy off and there is no evidence of the inertial waves that are apparent 
on large wires. The photographs of Siegel and Usiskin [87] and Siegel 
and Howell [85] show the same kind of process for greatly reduced 
gravity. The forces of surface tension and buoyancy remain important, 
however the process is slowed down and the effects of inertia are 
greatly reduced. Indeed, the soundtrack of a movie by Siegel and 
Keshock (associated with [88]) specifically notes the reduced inertial 
effects in boiling under low gravity. 
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Another question arises: If, at low gravity or for small wires, 

inertia becomes insignificant and the Zuber-Kutateladze mechanisms 
fail, what have the investigators who have given q max data for 
low R' actually observed? Without the inertial wave mechanisms it is 
hard to see how there could have been any q^ = ^ points to report. 

Ml O.X 

Experiment 

To answer the questions raised in the Introduction, Bakhru [20, 21] 
set out to measure the full q vs. &T curve for very small wires heating 
a variety of liquids. The apparatus used to do this is shown in Fig. 60. 
A small platinum wire, which served as both a resistance heater and a 
resistance thermometer, was suspended in the liquid of interest and 
boiling was observed at successive heat fluxes. The temperature of the 
wire was calculated from the resistance which in turn was computed 
from the ratio of voltage to current, using the method detailed by 
van Stralen and Sluyter [89]. Since complete details of the experimen- 
tal method are given by Bakhru [20], we shall only list a few major 
features of the tests here: 

The wires were cleaned in soap and then rinsed in the test liquid. 
Reagent grade liquids were used in all cases. During actual observa- 
tions the preheater, which was used to maintain saturation conditions, 
was momentarily switched off. Since the wires would melt during 
atmospheric pressure runs in water, the water runs were all made at 
pressures in the neighborhood of 75 mm Hg abs. The maximum pro- 
bable error in q was found to be 2-1/2%. The maximum probable 
error in At varied from 36% or about 3°F at the very lowest At's to 
1.8% or about 7°C at the highest At's. 

Bakhru's complete raw data are given in the form of 13 boiling 
curves in Figs. 61 through 72. These have been arranged in order of 
increasing R' from 0.0076 up to 0.0806. In each case data are pre- 
sented for both increasing and decreasing heat flux to expose any 
hysteresis effects. Figure 63 combines data for two wires under 
identical conditions to assure reproducibility of results. 

The major significance of Bakhru's study is readily apparent 
from these curves. As At is increased on the smaller wires, q rises 
monotonically without passing through a maximum and minimum. Only 
as R' increases to about 1/10 does the conventional boiling curve begin 
to re-establish itself. 
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Fig. 60 Schematic representation of experimental apparatus. 
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Fig. 61 Boiling curve for .0254 platinum wire in 
benzene. 
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Fig. 62 Boiling curve for .0254 platinum wire in 
acetone. 
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Fig. 63 Bolling curve for .0254 mm platinum wire 
in methanol. 


TTT 


1 1 — T 


TTT - 


i — i r 


.0254 mm pt wire in methanol 
saturated at I atm 
R' = . 00806 


.nucleation starts- 
jiucleation 


J_U 




film 
boilingj 

range of patchy_^ 
boiling from onset 
to full-blanketing 

st 


o increasing q | 


run 


• decreasing q 
□ increasing q 
■ decreasing q 

i i i i i 1 1 L 


► 


2 nc *run 


J i L 


144 
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Fig. 64. Boiling curve for .0254mm platinum wire in 
isopropanol. 
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Fig. 65 Boiling curve for .0508 platinum wire 
in water. 
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Fig. 66. Boiling curve for 076 mm platinum wire in 
water. 
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Fig. 67 Boiling curve for .0508 mm platinum wire 
in methanol. 
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Fig. 68 Boiling curve for .102 mm platinum wire in 
water. 
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Fig. 69 Boiling curve for .076 mm platinum wire in 
methanol. 
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Fig. 72 Boiling curve for .254mm platinum wire in 
methanol. 
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Still and high-speed motion pictures were made during the boiling 
process. For small R', the still pictures were similar to pictures 
presented by several previous authors for small wires at high heat 
fluxes ( e.g . [84] and [86]) as well as to those presented for large wires 
at low gravity ( e.g . [85] and [87]). The first nucleation resulted in a 
bubble which grew and spread horizontally until the wire was partially 
blanketed with a vapor patch. As q was increased, the fraction of the 
wire covered by vapor increased until the wire was totally blanketed. 

The photographs very strongly suggested that the blanketed portions 
of the wire were fully covered with vapor and that there was no liquid- 
to-heater contact within them. While these patches are properly to be 
called "film boiling" they exhibited none of the Taylor wave process. 
Instead, vapor removal was accomplished when buoyancy alone over- 
balanced surface tension, so that fairly large bubbles rose in an irre- 
gular way. 

The motion pictures provided details of patch growth and collapse. 
Figure 73 shows a tracing made from a 10-frame motion picture 
sequence during the growth of a patch on a wire in methanol. In the 
organic liquids that we used, patches would rapidly grow in this way 
and then collapse. Thus the vapor patches would flicker on and off 
from place to place on the wire. However, the patches were very 
stable in water. They would establish themselves and then spread 
steadily as the heat flux was increased, until the wire was totally 
blanketed. 

The boiling curves show that for R' < 0.01 the monotonic small- 
wire mechanisms are completely established. As R' is increased from 

0.01 to about 0.03, two phenomena appear: One is the onset of hydro- 

dynamic transition in the sense that nucleate boiling occurs first, and 
then undergoes a discontinuous transition to film boiling. The other 
phenomenon is a temperature overshoot on the rising heat flux curve. 

The resulting hysteresis in the curve becomes quite pronounced as R' 
is further increased to 0.08. The hydrodynamic mechanisms, of course, 
will only be fully re-established for R' values of about 0.12 or 0.15 — 
well beyond the range of the present tests. 

There are thus four distinct processes to consider in greater 
detail. They are: 

1. ) Natural convection from small wires. 

2. ) The inception of boiling 
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Fig. 73 A vapor patch propagation traced from a ten-frame motion 

picture sequence. Methanol on a 0.254 mm pt wire. R =.008, 
q =521,000 W/m 2 




3. ) The mixed-mode, or combined film-boiling and natural 

convection, regime. 

4. ) Film boiling from small wires. 

The nucleate and transition boiling regimes do not appear in the list 
since they do not occur when R' is small. 


Of these processes we shall concentrate on the last three. 

Bakhru's [21] consideration of the natural convection data presented in 
Figs. 61 through 72 revealed that they were generally consistent with 
the results of many investigators, as have been correlated by 
McAdams [SO], and that their scatter was generally comparable. A few 
points close to boiling inception fell somewhat above McAdams curve, 
probably owing to thermocapillary jetting from potential nucleation sites. 
This process is a Marangoni effect which has been discussed by a 
variety of investigators (see e.g. Trefethen [91], [92]). Such jetting 
was evident in our experiments and is disussed at greater length in 
[20]. For larger heaters, on which R' is small by virtue of low 
gravity, the conventional natural convection correlations for cylinders 
should apply perfectly well. 

The Inception of Boiling 

Double equilibrium of unstable nuclei. In 1961, Hsu [ 93 ] pointed 
out that when a nucleus bubble protrudes into liquid with a non-uniform 
temperature profile at a wall, there are two equilibrium positions that 
it might assume. Brown [ 94 ] subsequently described this process for a 
steady temperature profile, in the following way: 

The temperature profile at the wall is given in terms of the heat 
flux, q, and distance, y, from the wall, as 


T = T w - qy/k 


(119) 


The rate of change of temperature with pressure inside the bubble is 
given by the Clausius-Clapeyron equation as dT/dP = T v f /h . 

If we take the " r sa( . v^/h^ as nearly constant over the small^tempera- 

ture ranges of interest, integrate between T sa (. and T of interest, and 
finally introduce the Laplace equation (p - p = 2 a/r), we get the 
relation between the temperature in a bubbleand its radius, r, as 
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Now, with reference to Fig. 74 , let us consider a hemispherical 

nucleus protruding into a liquid whose temperature is given by equation 

(119) . At a point, r = y, for which the temperatures given by equa- 
tions (IIS) and (120) are equal, the cap of the bubble will be in equili- 
brium. However, the liquid around the equatorial region will be hotter 
than the gas inside, and evaporation will occur. This will cause the 
bubble to grow a little larger, forcing the cap into cooler liquid until 

a balance is struck between condensation at the cap and evaporation at 
the base of the bubble. Thus the equilibrium radius will be just a 
little in excess of the r for which the T's given by equations (119) and 

(120) are equal. 

Depending upon the parameters of equations (119) and (120), three 
situations can exist as suggested by Fig. 74. If there is no intersec- 
tion, the shape of any liquid-vapor interface that might exist in the 
cavity is determined, not by equations (119) and (120), but by such 
features of the cavity as the contact angle and the presence of non- 
condensible gases. If the curves simply fall tangent there is a single 
equilibrium radius, a little larger than r 2 . If the temperature distri- 
bution in the liquid near the wall is still steeper, there will then be 
two possible equilibrium points one of which will be a little larger than 
r 3 . The other equilibrium radius is smaller than the existing site in 
this case; hence it would have been triggered already in the cavity 
depicted. If the cavity were smaller, either site could exist. 

The double radii are stable with respect to slow changes of q 
(and of the temperature gradient) until the smaller radius becomes too 
small for the cavity. The larger radius will trigger if its size is some- 
how altered quickly enough that inertia will cause it to overshoot its 
equilibrium. Accordingly the group at M.I.T. named these larger 
nuclei, "metastable bubbles" and they noted experimental examples of 
their metastable behavior [94 , 95 ]. We shall have two occasions in the 
present discussion to point out their influence. 

12 

Nucleation hysteresis . In all cases, nucleation occured in 
Bakhru's tests at temperatures well in excess of that predicted by 
equations (119) and (120) for r^. Furthermore nucleation always began 
at temperatures a good deal larger than those at which it stopped. 


^We use the word "nucleation" to denote the inception of bubble growth 
whether it leads (as in the case of small R') directly to film boiling 
or (as in the case of larger R') to nucleate boiling. 
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Fig. 74. Equilibrium nucleus radii in a temperature 
gradient near a wal I. 
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Thus a "metastable bubble" undoubtedly gave rise to nucleation, and 
nucleation terminated when the temperature fell below the triggering 
temperature for the cavity. The disparity between the two points is 
called "nucleation hysteresis". 

Temperature overshoot hysteresis. The phenomenon of tempera- 
ture overshoot is well known. In prior studies, both Bakhru and Lien- 
hard noted [95, 96Q that when nucleation occurs the temperature 
sometimes drops back before continuing to rise with q. Then, the 
decreasing heat flux traverse usually drops to the right of the nuclea- 
tion heat flux and meets the convection regime with much less discon- 
tinuity. Several other authors, all the way back to Corty and Faust 
[97], have noted this hysteresis. In the present case the magnitude 
of the temperature overshoot diminishes with wire size as is shown 
in Fig. 75 — a plot of the temperature drop at inception, AT qv , against 
wire size. 

It is commonly argued (see, e .g . [97, 98]) that the overshoot occurs 
because the first nucleation triggers many adjacent sites and a sudden 
onset of widespread boiling cools the heater. Before the first boiling 
occurs, the larger nucleation sites which would normally have triggered 
inception, are initially filled with liquid and they do not function at low 
temperatures. But as the heat flux is decreased the sites, having filled 
with vapor and begun to function, will continue to function down to low 
temperatures, thus shifting the return traverse to the right and causing 
a hysteresis. That the temperature overshoot vanishes with decreasing 
wire size simply reflects the fact that there are fewer and fewer large 
sites on finely drawn wires. 

Thus the nucleation hysteresis should be the same for large wires 
at low gravity as for these small wires. But the temperature overshoot 
hysteresis, on the other hand, could not be expected to vanish on large 
wires at low gravity the way it does here. 

Mixed-Mode Boiling 

Vapor front propagation . Figure 73 suggests an important aspect 
in the growth of a vapor patch, namely that it spreads with a constant 
velocity. This behavior was reproducible in many motion picture runs 
for wires in methanol and it was evident in the other organic liquids, 
but not in water. Since this spreading behavior is a primary feature 
of the mixed-mode regime it should be treated in detail. 

Figure 76 shows the advancing vapor patch in a coordinate frame 
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Fig. 75 Influence of heater size on temperature 
overshoot. 
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Fig. 76 Steady-state representation of 
advance of a vapor patch. 



that moves with the speed of propagation, u. Thus the liquid-vapor 
interface is fixed while tine wire moves into it with a speed, u. Defi- 
ning 9 = T - T t , one can write the equilibrium wire temperatures 

under the patch and in the liquid as 0 and 0 , respectively. 

e b e 

The aim now is to determine the propagation or "triggering" 
temperature 0 O , that is consistent with an observed propagation speed. 
This will then help to reveal the mechanism of propagation. In parti- 
cular, it should reveal whether the interface propagates at the Leiden- 
frost temperature or at a nucleation temperature. 

The heat conduction equation for the wire is easily obtained from 
the familiar fin equation as long as the transverse Biot number, 
hR/k, is small: 

2 

u d8 d 0 , 2q 2h 0 

- — — = + — - 0 (121 

a dx .2 k R k R v 

dx w w 


where q here designates the local heat flux. It will be less in the un- 
blanketed portion of the wire since the temperature (and electrical resis- 
tance) are lower in that portion. Equation (121) will apply to either the 
blanketed portion of the wire (hereafter designated with the subscript, b) 
or to the unblanketed portion. The boundary conditions are 


e (°°)= e^. 
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b ©L 
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dx 
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9(0) = 0 b (O)=9 o (122) 


.x b =° 


The dependence of the propagation velocity upon system variables 
is given by u = u(a , q, q b , k^, R, h, h b , 0 Q ) where 0 0 and 0 are 

dependent variables since they equal q/h and q b /h b , respectively. There 
are thus 9 variables expressible in 4 dimensions, so the problem redu- 
ces to 5 independent dimensionless groups: 
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(123) 


The first group is a dimensionless triggering temperature. The 
second group is a kind of Peclet number involving the familiar fin 
constant, /2h/k w R. It compares the propagation capacity of the vapor 
front to the conduction capacity of the wire. The Biot number, Bi, is 
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typically on the order of 10 J in the present problem, and should there- 
fore exert no influence other than assuring us that radial conduction can 
be ignored in the total problem. Thus we anticipate a functional 
equation of the form: © 0 = 0 q (Pe, h b /h, q b /q). 


The analytical problem then reduces to a pair of coupled ordinary 
differential equations: 


0 " + p e 0’ - 0 = - 1 

0 " - p e /— 0 * - 0 = - h ^b _ _ 0 / 0 

b * h b b b h b q e b « 


(124) 


with boundary conditions: 


q ^ 

0(°°)= 1; © (») = JL. ; B'(0) = - ' -h 0' (0); 0(0)= © (0) = 9 

b qh V h b be 

b 


(125) 


where the dimensionless dependent variable is § = x</2h/k R or 

" x b v/ ^b. /k R. The solution is straight forward and the result is 

^ ^ D W 
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(126) 


Equation (126) will (by way of an example) now be compared with 
methanol data for a .0254 mm platinum wire (Figs. 61 and 73 ). To 
make the comparison, the blanketed and unblanketed conditions are taken 
to correspond approximately with the points at either end of the mixed- 
mode regime in Fig. 61 . Here, h = 4650 and h b = 2160 W/m 2 , so 
b^/b = 1/2.26. The temperatures at the end points are such that 
q b /q = 1.95. 

Figure 77 displays equation (126) for this q b /q and for three values 
b b /h: 0 , 1 / 8 , and the value of 1/2.26 given above. It is clear from 
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the curve that the relation between the dimensionless propagation 
velocity, Pe, and the triggering temperature, © , is not too sensitive 
to h b /h. In this case the average value of u was found to be about 
26. 2cm /sec (cf_. Fig. 15 where u = 25 cm/sec.) which gives Pe =2.48. 

The corresponding value of © G is 1.19, so the vapor front propagates 
by virtue of a triggering temperature only 19% higher than that which 
gives rise to the first nucleation. 

The triggering temperature for methanol is thus comparatively 
low — a typical nucleation temperature — and definitely not the 
Leidenfrost temperature. Figure 73 bears this out in that the front 
appears to be moving by triggering a succession of new bubbles imme- 
diately adjacent to it. The patch spreads until so much of the wire is 
insulated and heated-up that the resistance rises, throttling the current. 
This reduces the heat flux and the patch must again retreat. The 
result is the rapid, and apparently random, growth and collapse of 
patches on the wire. 

Non-propagation of the vapor front . One of the liquids used in 
Bakhru's study distilled water — did not show the propagation of the 
vapor patch. The patch (or patches) formed in water and spread in 
direct proportion to q. The flickering behavior was not present as it 
was in the various organic liquids. 

To see why this should be, one should recall (Figs. 61 through 72 ) 
that the nucleation hysteresis was far less pronounced in water than in 
the organic liquids. But nucleation hysteresis should depend on the 
presence of "metastable bubbles". Then, equating T as given by equa- 
tions (119) and ( 120 ), setting y = r and solving the result for r, one 
finds 


where 


hr = 1 
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The ratio of the two possible radii given by equation (127) will be real 
if nucleation occurs. Thus: 
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Equation (12S) is plotted in Fig. 78. Values of A based on the 
inception point for the smallest wire in each of the five liquids have 
been calculated and plotted on this curve. It turns out that the ratio 
r bi //p little ran 9 es between about 60 and 130 for the organic liquids 
but 9 it is only about 4 for water. Thus there are many large "metas- 
table bubbles” in the organic liquids. These tend to be triggered by 
the disturbances associated with the advancing vapor front, so the 
patch keeps over- running itself and then collapsing. Since any 
metastable bubbles in water differ relatively little from ordinary nuclei, 
the front can trigger only the nuclei very near the vapor— front and 
the advance is stable . 

Using Fig. 63 Bakhru found h b /h = 1/8 for water on a 0.0508 mm 
wire. Entering equation (11) with this value and a zero propagation 
speed ( i.e. , Pe = 0) he found ® 0 ~ 6.6 or (since 9 g ~ 21°C) 9 q ^ 240°C. 
The front can thus be at any temperature up to 9 = 240°C but it still is 
not triggering nuclei. Thus, without large metastable bubbles on the 
wire, only those nuclei immediately adjacent to the vapor-front will be 
triggered. Furthermore 9 Q is too low to trigger a Liedenfrost front. 

This is why the patch doesn't move in this case. 

Large cylinders at low gravity would behave like those described 
here with two possible changes: 1.) The Biot number might become 

large enough to make the heat conduction problem two-dimensional. 

This would require a considerably more complicated set of predictive 
relations than we have offered, possibly without great change in the 
results. 2.) With more nucleation sites on the cylinder, 9 q would 
probably tend to be lower than it is on small wires. 

Film Boiling 

Bromley [99] provided the first analytical prediction of heat transfer 
from cylinders in 1950. He found that: 

1/4 

N u /(Ra*) = 0.62 (130) 


whe re : 

Ra* = P (P - P ) h * g(2R) 3 /P k At (131) 

g f g fg 9 9 

In 1962, Breen and Westwater [59] provided a semi -empirical modifica 
tion of Bromley's equation. For 2R' <0.8 they recommended: 
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(132) 


Nu/(Ra*) 1/4 = 0.372(2R’) 1/4 + 0. 274/(2R’) 3/4 


Finally, in 1967 Baumeister and Hamill [58] derived the following 
relation 


= 0.373 j~(2R') + 3.68 + 1 ,083 . - 

!— J 


Bakhru plotted his film boiling data on Nu/(Ra*) vs . (2R') 
coordinates as shown in Fig. 79 , and included the preceding three 
expressions for comparison. Bromley's prediction is, of course, low 
for small R* and the other two expressions were developed with the 
objective of improving it. Equations (132) and (133) are a little high 
and a little low respectively, but either is in reasonable agreement 
with the available data. Since Baumeister and Hamill's expression in 
based upon an analytical rationalization that applies to wires of any 
size, the correlation in Fig. 79 should be valid for either large wires 
at low gravity or small wires under higher gravity. 

It is, of course, possible to fit an empirical curve to the data in 
Fig. 79 and obtain higher accuracy than either equation (132) or (133) 
will give. For example 

Nu/(Ra*) 1/4 ~ 0.7(2R') 1/3 (134) 

will fit the data nicely. 

Conclusions 

1. ) The local maxima and minima in the boiling curve which are 
dictated by hydrodynamic instabilities vanish for all R' < 0.01. The 
region 0.01 < R' < 0.15 is a transition regime in which the hydrodynamic 
mechanisms re-establish themselves. 

2. ) Nucleate boiling also vanishes when R' becomes small. 

3. ) Natural convection and film boiling on small wires (or on 
larger cylinders at low gravity) are predictable by conventional methods. 

4. ) The nature of a monotonic boiling curve between first nuclea- 
tion and full film boiling has been traced in detail and many of its 
features have been analyzed. This has been done for the constant 
electrical current case on small wires. The region is one in which 
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film boiling and natural convection coexist on the wire. If the para- 
meter A is on the order of a few hundredths or less, patches of 
film boiling flicker on and off over the wire. If A is larger than one 
or two tenths, the patches are steady and grow only with increasing 
heat flux. The analytical descriptions of the mixed-mode regime can 
generally be applied to large cylinders at low gravity as well. 

5.) Values of q max hitherto reported by many investigators 
(e.g. [84], [85], [86] for q max at R’ <0.03 have probably been either 
q at the inception of film boiling, or values of q for which the tempera— A 

ture first began to increase rapidly, or possibly values of q for which 
the heater first glowed red. In any case, such data are not meaning- 
ful since there is no longer a real local maximum to measure. 
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VIII. CONCLUSIONS 


Motivated by a need to know the influence of gravity on the boiling 
transitions, we have spent five years exploring the hydrodynamic mecha- 
nisms giving rise to these transitions. During the past three years 
these studies have been quite fruitful. In particular they have yielded 
major extensions of the hydrodynamic theory which, we believe, make 
a powerful problem-solving tool of the theory. The work has led to 
the following specific results: 

A. A broad scheme for predicting the peak pool boiling heat fluxes 
on a variety of heaters in low viscosity liquids . This material is 
developed in Chapters II and III. Predictions involving few or no 
empirical constants are developed for the following configurations: 

1. ) Horizontal flat plates without induced side flow: 

a) "infinite” plates (equation (15), Fig. 7) 

b) finite plates (equation (16), Figs. 8 and 9) 

2. ) Horizontal circular cylinders 

a) large (equation ( 31 ), Fig. 13) 

b) small (equation (37), Fig. 13) 

3. ) Spheres 

a) large (equation (38), Fig. 16) 

b) small (equation (30), Fig. 16) 

4. ) Horizontal ribbons, vertically oriented 

a) both sides uninsulated 

(i) large (equation (34), Fig. 18) 

(ii) small (equation (41), Fig. 18) 

b) one side insulated 

(i) large (equation (34), Fig. 19) 

(ii) small (equation (42), Fig. 19) 

In addition to making these predictions we have provided: 

5. ) Copious experimental verification of each prediction. 

6. ) Theoretical and experimental values of the vapor blanket 

thickness at the peak heat flux. (Fig. 12 and pg. 31). 
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7.) General guidelines for predicting on a variety of 

finite heater configurations. For these generalizations the 
reader should see pages 48 and 50 of the Conclusions of 
Chapter III. 

B . Theoretical descriptions of the, Taylor and Helmholtz stability 
mechanisms underlying the hydrodynamic theory of the boiling transi- 
tions . Solutions of these problems, although developed for application 
to boiling heat transfer, have potential application in any of a large 
variety of two-phase dynamical configurations. The two problems 
are: 

1. ) the prediction of Taylor instability in the top of the inter- 

face between a horizontal gaseous cylinder and the surroun- 
ding viscous liquid. The horizontal plane interface between 
a gas (below) and a viscous liquid (above) is obtained as 
the special case for large radius. The analytical solution 
is given in Figs. 23 through 26. The results show that 
viscosity slows the growth rate of Taylor waves and con- 
siderably increases their length. 

2. ) the prediction of the maximum velocity of gas for which a 

vertical jet in a viscous liquid remains stable. This result 
is given by equation (105) in which the three decimal-con- 
stants must be viewed as arbitrary. 

C. Experimental studies of film boiling. A variety of new data have 
made possible: 

1. ) a verification of the predicted Taylor wavelengths and fre- 

quences for viscous liquids. (Figs. 35 through 38). 

2. ) an empirical correlation of vapor blanket thickness during 

film boiling. (Fig. 28). 

3. ) an empirical correlation of an observed stretching of the 

wavelength at high vapor volume fluxes, independent of the 
influence of viscosity. (Fig. 45). 

4. ) verification of the lack of correlation between bubble fre- 

quency and wave "frequency" in film boiling, (pgs. 91 and 
97, and Fig. 44). 
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D. Peak and minimum pool boiling heat fluxes in viscous liquids. 

1. ) Item B-2, above, was used in the formulation of a theoreti- 

cal expression for q max in viscous liquids, (equation (105)) 

2. ) Original data provided verification of the viscous q,^^ 

theory. (See Figs. 50 and 51 and Conclusions on page 117). 

3. ) It is shown, analytically, that viscosity influences q m ^ n in 

complicated and conflicting ways (pg. 101). This problem 

is left unresolved. 

ET. Hydrodynamic theory of high current-density electrolysis . Like 
the pool boiling heat flux, the current density during electrolysis 
passes through a peak, after which a film electrolysis region is 
established. It is predicted, by Bhattacharya, and verified experimen- 
tally, that this peak corresponds with Moissis and Berenson's "first 
hydrodynamic transition". A prediction for the peak current density 
(equation (117))is formulated and verified experimentally (Fig. 58). 
More detailed conclusions are given on pg. 137. 

F . The vanishing of the hydrodynamic mechanisms at low gravity 
and for small sizes of heaters. It is shown by the extensive experi- 
ments of Bakhru that the hydrodynamic stability mechanisms vanish 
as R' is reduced from 0.15 to about 0 = 01, For R’ < 0.01 the con- 
ventional boiling curve with its well-known maxima and minima is 
replaced with a monotonic curve. This curve passes from natural 
convection directly into film boiling without a nucleate or transitional 
boiling regime. Full conclusions relating to the prediction of the 
heat flux in this case are given in the conclusions on pages 168 and 
170. 


G. "Spinoff". A variety of studies completed under the present 
grant, but not immediately related to the hydrodynamic theory, are 
excluded from this report. These studies of condensation, convection, 
fins, and computer coding are described briefly on pages 2 and 3, and 
full references are given to published work. 

The initial aim of determining the influence of gravity has been 
fully met in this work since every one of these results includes 
gravity in the analysis or in the non-dimensional correlation of data. 
The influence of gravity is thus made clear in all cases. 
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APPENDIX A NOMENCLATURE 


A 


h 


area of heater (or of cathode in the context of elec- 
trolysis) 


A 

J 

B 

Bi 

Bo 



cross-sectional area of vapor jets escaping from A^ 

ratio of viscous to surface tension force in vapor 
blanket, qM g /(P g h f * a) 

Biot number, hR/kyy in context of Chapter III and 
hR/2k w in context of Chapter VII. 

the Bond number, ratio of buoyant to surface tension 
forces, R§(P f - P ) g/c = R’ c 2 

constant in expression, equation (111) 

constant in v expression, equation (112) 

min . 
cyl 

fraction of axial distance, on cathode or heating 
element, occupied by gas or vapor jets 


c 


speed of a disturbance in a liquid-vapor interface 


c (T) specific heat of heater material, a function of T 

c y specific heat of vapor 


D 


b 


bubble diameter at departure from cathode or heating 
element 



E 

E 

c— b 


f( ) 


g 
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depth of liquid and vapor phases, respectively 
voltage drop between cathode and anode 
voltage drop from cathode to nearby liquid bulk 
any function of 

number of bubbles released per second 

actual gravity (or body) force acting on a heater or 
cathode 


earth -normal gravity 


vertical dimension of a horizontal ribbon 
latent heat of vaporization 

latent heat of vaporization corrected for sensible heat. 



heat transfer coefficient 


electrical current 

f=T 

dimensionless wave number, k ^/[g(Pf - P ) ]. Also 
denotes dimensionless group defined in Fig. 28 

thermal conductivity of liquid, vapor, and heater, 
respectively 

a characteristic dimension (= H or R in certain of the 
present applications). Also used in Chapter VI to 
denote spacing between adjacent vapor jets 

liquid viscosity parameter, defined in equation (75). 

Also used to designate molecular weight (of Hg) in 
Chapter VI 

Borishanski number, defined in equation (11) 
number of vapor jets on a heater of area, A 

h 

Nusselt number, 2hR/k 

number of electrons (equal to 2) per Hg atom formed by 
electrolysis 

nucleation site density 

absolute pressure (at surface of heater in context of 
Chapters IV and V). Also used in context of Chapter 
III to denote length of perimeter of the cross-section 
of a long slender heater 


dimensionless velocity, or Peclet number, defined by 
equation (123) 


Pg pressure in a gas jet 

P f > P g perturbation pressure in liquid and gas, respectively 

P Q total perturbation pressure 


q average heat flux on an entire heater. Also heat flux 

on the unblanketed portion of the heater when used 
with q b in Chapter VI 

^ peak heat flux for a particular heater geometry and 

F for an "infinite" flat plate, respectively 

a characteristic heat flux defined by equation (1), 

2 equal to Zuber's prediction for infinite horizontal flat 

plates 


^min* ^min minimum heat flux on a cylindrical heater and on an 
F infinite flat plate, respectively 


q p radiant heat flux 

R radius of spherical or cylindrical heating element or 

of cylindrical cathode 

Ra * modified Rayleigh number defined in equation (131) 

R c radius of a cylindrical heater corrected for the 

vapor blanket thickness around it 


P® Reynold's number as defined below equation (94) in 

the context of Chapter V, and as defined in equation (86) 
in the context of Chapter IV. 

Rj radius of an escaping vapor jet 

r radius of a bubble 


T tern pe ratu re 

temperature overshoot, T at nucleation minus minimum 
T for q beyond nucleation 
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sat 


w 


t 

u f 

u t 

u 


m 


U 


’m 


U 


H 


u 


U g* u f 


saturation temperature of a liquid 
temperature of heater wall 
time 

primary velocity of the liquid in a liquid column 

maximum primary downward velocity in the liquid 

primary velocity of the gas in a gas jet 

maximum primary velocity of the gas 

vapor velocity in a jet, for which jet becomes 
Helmholtz unstable 

speed of vapor patch propagation. Also used to denote 
perturbation velocity of gas in x-direction in Chapter V 

perturbation velocity of gas and liquid, respectively, 
in the x-direction. Used in Chapter IV 


V dimensionless gas viscosity papameter, defined in 

equation (106) 

v perturbation velocity of gas in the y-direction. Used 

in Chapter V 


v vapor or gas volume flux 

v , v perturbation velocity of liquid and gas, respectively, 

f 9 in the y-direction. Used in Chapter IV 

v , v , v volume flux at the "second" (or Helmholtz unstable) 

max min sc 

transition, at the minimum film electrolysis or boiling 
transition, and at the "first" (or isolated-bubbles to 
slugs-and-columns) transition, respectively 

x, x axial distance coordinate; subscript b denotes reversed 

b coordinate under the blanketed portion of the wire 


distance measured normal to axial direction 
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Greek Letters 


a 

0 


6 

A 


'T 


At, Ar 


thermal diffusivity of heater 


contact angle or angle between liquid-vapor (or gas) 
interface and heater (or cathode) surface, measured on 
the liquid side 


dimensionless density, (p - p ) / (p -f p ) 

f 9 f 9 

vapor blanket thickness near the peak heat flux 


dimensionless vapor blanket thickness in film 
d /R. Also designates 6/g(p^ - p )/o 

g g 


boiling. 


temperature difference between a heater wall and a 
liquid 


increments of time and sphere radius in numerical 
computation 


n 


0 


9 



e 

G 

e 

0 

o 


A 


d 


X 


X 


c* 


ordinate of an interface measured from its mean 
position 

dimensionless tern peratures, 6/0 and 9 /0 , 
respectively e ° e 


temperature measured above T 

sat 

equilibrium temperature in the wire, q/h 
vapor patch triggering temperature 


dimensionless wavelength, (X/2rr)/g(p 

Also denotes a dimensionless function, 
in Chapter VII 


- P ) / 3C = 1/K /~ 3 ~ , 

2chv T /ATh„ k 
rg sat rg 


dimensionless "most susceptible" wavelength 
wavelength 

the "critical" or minimum and "most susceptible" or 
most rapidly growing wavelengths, respectively 
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Helmholtz unstable wavelength 

viscosity of liquid and vapor, respectively 

dimensionless axial parameters defined after 
equation (125) 

density of saturated liquid or of electrolyte 

density of saturated vapor or of gas generated at 
cathode 

density of heater material 

density of gaseous hydrogen (special case for p in the 
context of electrolysis) ^ 

surface tension between a liquid and its vapor 

perturbation shear stress in liquid and vapor, 
respectively 

perturbation potential function as defined in equation (57) 

perturbation function as defined in equation (52) 

dimensionless frequency, uu /a/g 3 (p - p ) 

dimensionless "most susceptible” frequency 

angular frequency (complex). Denotes ~^.(W) in Chapter IV. 

maximum or "most susceptible" angular frequency 

"most susceptible" angular frequency on a flat plate 
when both the liquid and gas are considered inviscid, 

4/ 

0.620 g 3 (Pjr — P g )/C 
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General Subscripts 


b denoting conditions related to the portion of heater 

blanketed by vapor 

c denotes a critical, or maximum stable, condition 

cyl denotes horizontal cylinder 

F "infinite'* flat plate 

General Superscript 

denotes a length multiplied by /g(p^ - p )/a 
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